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- ABSTRACT?-/ The  design  of  vehicles  with  particular  operating  character¬ 
istics  on  rough  terrain  requires  accurate  descriptors  of  terrain  roughness. 

At  present  a  one  parameter  description  of  surface  roughness  is  available 
that  uses  the  standard  deviation  of  a  detrended  elevation  profile.  In  order 
to  introduce  the  frequency  content  of  the  terrain  into  the  roughness  de¬ 
scription,  a  three  parameter  model  of  the  surface  roughness  power  spectrum 
is  developed  whose  parameters  are  determined  from  the  values  of  the  standard 
deviation  of  displacement,  slope  and  curvature  obtained  from  a  detrended 
elevation  profile.  The  three  parameter  model  predicts  five  distinct  types 
of  power  spectra  which  can  be  used  to  classify  terrain  roughness  including 
the  cases  where  periodic  features  are  present.  For  some  types  of  terrain 
it  is  not  necessary  to  use  detrended  data,  and  for  these  cases  the  power 
spectra  are  developed  using  an  undetrended  formalism. 

The  three  parameter  roughness  model  can  be  used  to  predict  the  power 
absorbed  by  the  driver.  The  effects  of  vehicle  geometry  including  wheel 
size,  wheel  spacing  and  track  length  are  introduced  in  a  natural  way  through 
the  roughness  model.  Vehicle  dynamics  effects  including  the  location  of  the 
driver  relative  to  the  center  of  mass  of  the  vehicle  are  introduced  by  trans- 
mission  ^unctions .  A  comparison  with  experimental  absorbed  power  data  is 
made .  <r — 


PART  I:  INTRODUCTION.  The  dynamic  response  of  a  specified  type  of  wheeled 
vehicle  operating  on  a  terrain  depends  on  the  vehicle  speed  and  the  terrain 
roughness.  For  satisfactory  operation  of  both  vehicle  and  driver,  the  dynamic 
response  of  the  vehicle  must  be  predictable  for  a  given  type  of  terrain.  To 
achieve  this  prediction  capability,  the  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES)  was  requested  to  develop  improved  methods  of  describing 
terrain  and  use  them  to  calculate  the  resulting  vehicle  and  driver  motion. 


A  description  of  terrain  roughness  is  especially  important  for  the  de¬ 
velopment  of  new  vehicles  because  ride  tests  on  various  types  of  terrain  are 
required  to  evaluate  new  design  concepts.  An  accurate  surface  roughness 
description  is  required  to  separate  the  effects  of  a  new  vehicle  design  from 
the  effects  introduced  by  the  surface  roughness  properties  of  a  test  site. 
Both  the  frequency  content  and.  the  amplitude  of  an  elevation  profile  will 
influence  the  ride  quality  of  a  vehicle,  and  an  effect  observed  during  a 
vehicle  ride  test  may  be  due  to  a  peculiarity  of  the  terrain  roughness  rather 
than  a  specific  vehicle  effect. 

The  development  of  new  vehicle  designs  is  helped  by  a  thorough  under¬ 
standing  of  how  changes  in  the  geometry  and  dynamics  of  a  vehicle  will  inter¬ 
act  with  the  terrain  roughness  to  produce  changes  in  the  ride  quality. 


Especially  important  are  the  effects  of  geometry  including  wheel  size,  wheel 
spacing,  driver  location  and  track  length,  and  the  effects  of  dynamical  param¬ 
eters  including  spring  constants,  damping  constants  and  the  weight.1  The 
effects  of  soil  conditions  are  not  treated  although  they  may  be  important  to 
the  vehicle  response. 

« 

Description  of  Terrain  Roughness 

A  description  of  surface  roughness  begins  with  a  measurement  of  an  ele¬ 
vation  profile,  Figure  la.  Elevations  are  measured  at  small  intervals 
(usually  one  foot  or  one  half  foot)  over  a  distance  of  several  hundred  feet. 
Numerical  differentiation  gives  the  slope  and  curvature  profiles.  The  mea¬ 
sured  elevations  include  gradual  changes  of  elevation  over  a  long  distance. 

In  some  cases  this  nonstationary  character  of  the  measured  elevation  profile 
may  have  to  be  removed  by  a  detrending  procedure  in  order  for  the  data  to  be 
processed.  In  other  cases  the  measured  elevation  data  can  be  used  directly. 

Roughness  can  be  described  as  a  stationary  random  process  if  the  statis¬ 
tical  properties  of  the  profile  height  do  not  change  with  position  along  the 
profile.  To  remove  nonstationary  trends  from  the  data  a  detrending  procedure* 
is  applied  to  the  data  which  removes  the  very  long  wavelengths  (the  trend), 
see  Figure  lb.  The  mathematical  procedure  for  obtaining  detrended  data  is 
given  by2 

00 

.  Fd(x)  =  F(x)  -  jx  /  [F(x  +  a)  +  F(x  -  a)]e"aM  da  (1) 

o 

where  F,(x)  =  detrended  function 
a 

F(x)  =  original  function 

\  =  filter  constant  (detrending  constant) 

Computer  program  RRFN  has  been  developed  to  calculate  the  detrended  data  from 
the  original  input  data.1 

After  removing  the  trend  the  question  still  remains  as  to  the  frequency 
constant  of  the  terrain  roughness.  This  is  important  because  a  dynamical 
system  such  as  the  spring  support  of  a  vehicle  has  a  resonance  frequency,  and 
this  resonance  may  be  excited  if  the  vehicle  runs  over  a  given  terrain  at  a 
specified  speed.  Therefore  it  is  of  importance  to  know  the  wavelength  content 
of  the  terrain  roughness  because  this  will  determine  the  frequency  content  of 
the  input  displacement,  velocity,  and  acceleration  on  the  wheels  of  the  moving 
vehicle.  The  wavelength  content  of  a  terrain  displacement  profile  can  be 
described  by  a  power  spectral  density  function. 

The  power  spectral  density  (PSD)  gives  the  spectral  density  of  the 
variance  of  a  random  process.2'4  Therefore, 


ct2  =  V  =  /“  PSD  d« 
o 


(2) 


*  "Detrend"  means  to  remove  the  trend  of  the  elevation  variations. 
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where  a  =  RMS  value  of  the  random  process  =  standard  deviation  for  zero 
average  value 
V  =  variance 

Q  =  spatial  frequency  =  reciprocal  wavelength  =  number  of  cycles /me ter 
In  this  paper  the  RMS  value  and  the  standard  deviation  will  be  used  inter¬ 
changeably  because  the  baseline  can  be  chosen  so  that  all  average  values  are 
zero.  The  PSD  gives  a  measure  of  the  frequency  content  of  the  random  process. 
For  example,  the  surface  displacement  profile  PSD  is  a  measure  of  the  wave¬ 
length  content  of  this  profile. 

One  Parameter  Model 

A  descriptor  commonly  used  in  terrain  roughness  models  is  the  standard 
deviation  (RMS)  of  displacement  obtained  from  a  measured  elevation  profile. 

The  one  parameter  roughness  power  spectrum  can  be  calculated  from  the  value 
of  the  RMS  displacement.  The  one  parameter  power  spectra  of  the  surface  dis¬ 
placement  and  slope  can  be  written  as2*5 

P.(Q)  =  —  (3) 

d  Q2 

V«)  =  »^2n2[s-"reSM,)]2  V°>  “> 

where  P^  =  PSD  of  displacement  profile 
P^  =  PSD  of  slope  profile 
C  =  roughness  parameter 

AL  =  interval  of  elevation  of  profile  measurement 
For  the  undetrended  case  the  RMS  displacement  integral  (2)  does  not  converge 
with  P^  given  by  equation  (3),  however  the  RMS  slope  integral  (2)  does  converge 
with  P°  given  by  equation  (4).  For  the  single  parameter  model  a  theoretical 
expresiion  for  the  RMS  value  of  the  displacement  can  only  be  determined  for 
detrended  data,  and  the  theoretical  need  for  detrending  is  that  it  allows  the 
elevation  power  spectrum  integral  to  be  evaluated.  For  this  reason  the  value 
of  the  parameter  C  must  be  determined  from  detrended  data.  The  spectral  win¬ 
dow  function  occurring  in  equation  (4)  introduces  the  interval  of  measurement 
of  the  elevation  profile. 

The  power  spectrum  for  a  detrended  elevation  profile  can  be  written  as2 


Pd,DET(fi,A) 


- (2ttA)4Q4 — ^  p  (Q) 

|1  +  (2nMl)2] 

=  <r2Pd(«) 


(5) 


where  the  detrending  function  <)>  is  defined  by 


$  =  1  +  (2nAfi)'2 


(6) 
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Since  the  RMS  integral  of  the  slope  power  spectrum  converges,  Van  Deusen  uses 
the  undetrended  slope  power  spectrum  given  by  equation  (4)  to  calculate  the 
RMS  value  of  the  slope.  The  integrals  presented  are2 


o2  (Undetrended) 
d 


a2  (Detrended ,A) 


a2  (Undetrended) 
s 


J  P,(0)dft  ‘ 
o  d 


Diverges 


f  Pd,DET(0,X)dfi  =  2 


rPs(Q)dO  = 


(7) 

(8) 

(9) 


where  =  standard  deviation  of  terrain  displacement 
cr  =  standard  deviation  of  terrain  slope 

A  determination  of  the  RMS  displacement  for  detrended  data  (by  computer 
program  RRFN,  for  example)  immediately  determines  the  constant  C,  and  there¬ 
fore  the  power  spectrum  for  the  undetrended  data  is  calculated.  Two  basic 
theoretical  concepts  are  introduced  in  the  one  parameter  model:  (a)  the 
detrending  function  which  allows  the  elevation  RMS  value  to  be  calculated 
analytically,  and  (b)  the  slope  spectral  window  function  which  introduces  the 
measurement  interval  of  the  elevation  profile.  The  value  of  a  can  be  used  to 
calculate  C  directly  from  undetrended  data  using  equation  (9). 

Three  Parameter  Model 

A  comparison  of  the  one  parameter  monotonically  decreasing  power  spectrum 
Pj  =  C/ft2  with  the  power  spectrum  values  obtained  by  calculating  the  Fourier 
transform  of  the  autocorrection  function  of  an  elevation  profile  suggests  that 
the  one  parameter  representation  is  valid  only  for  high  spatial  frequencies 
of  the  terrain  roughness.  Also,  quasi-random  terrain  areas  such  as  plowed  or 
disked  agricultural  fields  cannot  be  described  by  a  one  parameter  model  be¬ 
cause  a  pronounced  peak  occurs  in  the  measured  power  spectrum  at  a  frequency 
corresponding  to  the  width  of  the  plowed  furrows. 

A  more  complete  description  of  surface  roughness  is  needed,  one  that  goes 
a  step  beyond  a  one  parameter  characterization  that  is  based  on  a  single  RMS 
value  of  surface  displacement.  This  paper  introduces  a  three  parameter  rough¬ 
ness  power  spectrum  model  that  can  describe  terrain  areas  having  periodicities, 
and  gives  a  better  description  of  the  roughness  power  spectrum  in  the  long 
wavelength  region. 

The  three  power  spectrum  parameters  are  determined  from  the  values  of  the 
standard  deviations  of  the  displacement,  slope  and  curvature  of  a  measured 
elevation  profile  (Figure  la).  The  power  spectrum  determined  by  these  param¬ 
eters  is  found  to  be  of  five  distinct  types  of  which  the  Van  Deusen  form  Cft"2 
is  one.  The  displacement  power  spectrum  P^ft)  is  found  to  be  a  complicated 
function  of  frequency  that  includes  spectral  windows  of  definite  bandwidth. 
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The  resulting  frequency  channels  of  power  transmission  and  power  suppression 
may  be  responsible  for  the  anomalous  dynamic  response  of  some  vehicles  on  a 
given  terrain. 

Actual  terrain  locations  where  military  trucks  and  track  vehicles  are 
expected  to  operate  can  vary  considerably  in  quality  (frequency)  as  well  as 
their  RMS  values  of  elevation.  Terrain  types  such  as  trails,  cultivated 
fields,  virgin  terrain,  test  courses,  etc.,  are  expected  to  be  encountered  by 
military  off-road  vehicles.  An  attempt  is  made  to  classify  terrain  roughness 
in  terms  of  the  five  basic  types  of  roughness  power  spectra  by  associating 
each  terrain  site  with  one  of  the  five  types  of  power  spectra.  The  type  of 
power  spectrum  associated  with  a  terrain  site  will  determine  the  nature  of 
the  dynamic  response  of  a  vehicle  moving  over  that  terrain.  In  this  way  the 
expected  dynamic  response  of  a  vehicle  can  be  related  to  a  specific  terrain 
site.  A  computer  program  TERR  was  developed  to  accomplish  this  work.1 

Absorbed  Power 

Driver  fatigue  can  be  related  to  the  power  absorbed  by  a  driver  during 
the  operation  of  a  vehicle.  Absorbed  power  is  a  physiological  concept  and  by 
definition  it  is  related  to  the  acceleration  power  spectrum  measured  at  the 
drivers  seat.6  Experience  has  shown  that  a  driver  can  absorb  no  more  than 
six  watts  of  power  and  be  physically  able  to  drive  a  vehicle  for  an  extended 
period  of  time.7  It  is  of  military  interest  to  be  able  to  predict  the  vehicle 
speed  that  corresponds  to  the  six  watt  absorbed  power  level,  and  to  determine 
the  terrain  roughness  descriptors  required  to  do  this. 

At  present  the  six  watt  speed  is  estimated  from  the  RMS  value  of  surface 
displacement  as  determined  from  a  detrended  elevation  profile.8-13  This 
single  parameter  terrain  roughness  description  is  sometimes  found  to  be  inad¬ 
equate  because  the  frequency  content  of  the  ground  surface  elevation  is  not 
considered.  The  three  parameter  terrain  roughness  description  was  developed 
in  part  to  improve  the  method  of  determining  the  six  watt  speed  by  including 
the  effects  of  frequency  content. 
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PART  II:  THEORETICAL  ROUGHNESS  MODELS 


INTRODUCTORY  REMARKS.  Terrain  areas  can  be  natural  or  man-made  or  a  com¬ 
bination  of  both  types.  Natural  terrain  elevations  appear  to  be  random,  but 
some  features  such  as  windblown  sand  dunes  have  an  apparent  regularity.  Man¬ 
made  areas  such  as  roads,  trails,  plowed  and  disked  fields,  etc.,  can  contain 
apparent  regularities  although  their  elevation  profiles  are  still  character¬ 
ized  as  random  data.  Plowed  fields,  for  example,  exhibit  a  quasi-sinusoidal 
pattern  and  will  have  a  power  spectrum  that  peaks  at  a  spatial  frequency  that 
corresponds  to  the  quasi-sinusoidal  pattern.  This  case  clearly  cannot  be 
described  by  a  one  parameter  power  spectrum  of  the  form  CQ'2  which  is  a 
monotonically  decreasing  function  of  spatial  frequency.  This  part  of  the 
paper  develops  a  three  parameter  power  spectrum  model  that  can  be  used  to 
describe  random  terrain  roughness  elevation  which  include  regularities,  man¬ 
made  or  natural. 

The  power  spectrum  coefficients  of  the  three  parameter  model  are  calcu¬ 
lated  from  the  RMS  values  of  elevation,  slope  and  curvature  that  are  obtained 
from  a  detrended  elevation  profile  or  from  an  undetrended  elevation  profile 
when  possible.  The  values  of  these  three  parameters  can  be  used  to  classify 
terrain  roughness. 

FORM  OF  THE  THREE  PARAMETER  POWER  SPECTRUM.  This  section  suggests  an 
appropriate  form  for  the  three  parameter  roughness  power  spectrum.  The  choice 
of  the  form  of  the  three  parameter  power  spectrum  model  is  a  natural  extension 
of  the  Van  Deusen  one  parameter  model  CQ"2 .  Figure  2a  shows  that  the  PSD  of 
the  undetrended  and  detrended  displacement  profile  is  described  by  the  form 
CQ'2  only  in  the  region  of  large  Q  (short  wavelengths).  For  small  Q  (long 
wavelength)  the  PSD  deviates  from  that  predicted  by  CQ'2.  Additional  param¬ 
eters  are  required  to  describe  the  deviation  from  CQ'2  in  the  region  of  small 


that  are  added  to  the  CQ'2  term.  In  this  way  the  modification  of  the  PSD 
becomes  important  only  for  small  Q  and  approaches  the  form  CQ'2  for  large  Q. 

The  PSD  for  undetrended  surface  displacement  can  be  written  as  a  general¬ 
ization  of  CQ'2  in  the  form  of  a  power  series  as  follows 


PJ  =  CQ'2  +  DQ"3  +  EQ'4  +  ••• 

a 


(10) 


For  practical  calculation,  the  power  series  must  be  cut  off  at  some  finite 
term.  If  two  parameters  are  chosen  to  describe  in  the  form 

Pd  =  CQ'2  +  DO'3  (11) 

then  the  measurement  of  both  a.  and  a  is  required  to  determine  C  and  D.  If 

d  s 

three  parameters  are  chosen  to  describe  P^  in  the  form 
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P.  =  Cft-2  +  Dft*3  +  Eft'4 

a 


(12) 


then  CT^,  Og  and  ac  would  be  required  to  determine  C,  D,  and  E,  where 

cr =  standard  deviation  of  curvature  (second  derivative  of  elevation  with 

regard  to  space  coordinate).  In  this  paper  the  term  curvature  will  be  used 
to  describe  the  second  derivative. 

Polynomial  models  beginning  with  the  term  Aft”1  such  as 

P.  =  AT)'1  +  Cfi-2  +  DQ-3  (13) 

d 


are  not  valid  for  two  reasons:  (1)  they  do  not  approach  the  Van  Oeusen  form 

Cfi-2  for  large  ft,  and  (2)  the  integral  f° P.<|>“2dft  for  detrended  data  has  a 

o  d 

logarithmic  singularity  and  cannot  be  evaluated. 

The  short  wavelength  (high  frequency)  part  of  the  power  spectrum  is  de¬ 
scribed  by  the  parameter  C,  while  the  long  wavelength  (low  frequency)  part  of 
the  powei  spectrum  is  described  by  the  parameter  E.  The  middle  ranges  of 
frequency  are  determined  by  the  parameter  D.  The  values  of  the  power  spectrum 
parameters  C,  D  and  E  are  determined  by  the  type  of  terrain  elevation  profile 
that  is  measured,  and  these  parameters  can  be  used  as  descriptive  numbers  to 
classify  terrain  roughness. 

CLASSIFICATION  OF  TERRAIN  ROUGHNESS.  This  section  shows  how  a  terrain 
roughness  classification  can  be  developed  in  a  natural  way  from  the  algebraic 
signs  of  the  parameters  C,  D  and  E.  Terrain  roughness  can  be  classified 
according  to  the  shape  of  the  power  spectrum  associated  with  the  natural  ele¬ 
vation  profile.  The  shape  of  a  theoretically  predicted  roughness  power  spec¬ 
trum  is  determined  by  the  values  of  the  parameters  C,  D  and  E.  These  param¬ 
eters  will  be  determined  from  the  values  of  a^,  Os  and  ct^  associated  with  the 
elevation  profile. 

A  priori,  the  values  of  the  three  parameters  can  be  positive  or  negative 
yielding  a  total  of  2  «  2  x  2  =  8  possible  types  of  power  spectra.  But  the 
case  C  <  0,  D  <  0  and  E  <  0  must  obviously  be  excluded  because  the  power  spec¬ 
trum  must  be  positive  at  least  in  some  frequency  range.  This  leaves  seven 
possibilities.  But  of  these  seven  cases  two  are  redundant  and  only  five 
physically  distinct  types  of  elevation  power  spectra  are  possible  for  the 
three  parameter  model.  These  are  shown  in  Figure  2b.  The  basic  forms  of 
Type  2  and  Type  3  remain  unchanged  for  D  >  0  or  D  <  0  so  that  either  sign 
still  represents  only  one  spectral  type,  and  this  redundancy  yields  only  five 
distinct  spectral  classes. 

The  one  parameter  power  spectrum  Cft'2  is  seen  to  be  a  special  case  of  a 
Type  1  power  spectrum  with  C  >  0,  D  =  0  and  E  =  0.  Type  1  power  spectrum  con¬ 
tains  all  frequencies  with  little  power  at  the  high  frequencies.  Type  2  has 
no  power  at  high  frequencies.  Types  3  and  A  have  pronounced  peaks  at  definite 
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frequencies  and  can  be  used  to  describe  a  large  scale  rolling  type  of  terrain. 
Type  4  with  D2  <  4CE  gives  a  power  spectrum  which  is  all  negative  and  so  this 
case  must  be  excluded.  Type  5B  with  D2  >  ACE  describes  the  situation  where  a 
high  frequency  periodicity  occurs  within  the  terrain,  such  as  in  the  case  of 
plowed  or  disked  fields.  Type  5A  with  D2  <  4CE  can  be  considered  to  be  a 
special  case  of  Type  1  where  the  power  is  positive  for  all  frequencies. 

a 

Power  spectra,  by  definition  as  the  Fourier  transform  of  an  autocorrela¬ 
tion,  are  always  positive.  Therefore  only  the  frequency  ranges  where  P^fQ) 
is  positive  in  Figure  2b  are  physically  acceptable.  The  regions  where 
P^(Q)  become  negative  must  be  regarded  as  physically  unacceptable  for  power 
transmission  and  P^(Q)  is  set  equal  to  zero  here.  In  this  manner  it  is  seen 
that  the  roughness  power  spectrum  P^O)  is  not  in  general  a  simple  function 
but  rather  may  contain  cutoff  frequencies  (fi0  and  Qj)  which  produce  "windows" 
(bright  areas)  and  zero  power  areas  (dark  regions)  of  definite  bandwidths. 
Therefore  only  in  distinct  frequency  regions  may  power  be  transferred  to  a 
moving  vehicle.  The  type  of  power  spectrum  exhibited  by  a  given  site  is  de¬ 
termined  from  measured  values  of  a^,  Og  and  Oc  from  which  the  values  of  C,  D 
and  E  are  obtained. 

The  power  spectra  used  are  one-sided  because  they  are  defined  only  for 
positive  frequencies.  Sometimes  two-sided  power  spectra  are  used  which  are 
defined  for  positive  and  negative  frequencies.  The  two-sided  power  spectra 
can  be  obtained  as  mirror  images  of  the  one-sided  power  spectra.  Power 
spectra  must  be  symmetrical  about  the  zero  frequency  axis  because  they  are 
defined  as  Fourier  transforms  of  autocorrelation  functions  which  themselves 
are  symmetrical.  The  power  spectrum  model  given  in  (12)  can  be  used  for  nega¬ 
tive  Q  by  choosing  | Q| 3 ,  absolute  value  of  0,  in  the  regions  of  negative  Q. 

The  type  of  power  spectrum  obtained  from  an  elevation  profile  will  depend 
on  the  length  of  the  measured  profile.  Therefore  the  classification  of  ter¬ 
rain  areas  will  be  scale  dependent,  and  a  comparison  of  two  terrain  areas 
should  be  made  on  the  basis  of  equal  elevation  profile  lengths.  This  paper 
does  not  investigate  scale  dependence. 

The  five  possible  types  of  roughness  power  spectra  are  associated  with 
the  choice  of  a  three  parameter  power  spectrum  model.  The  choice  of  a  greater 
number  of  parameters  would  produce  more  complicated  power  spectrum  shapes. 

For  example  a  four  parameter  model  would  exhibit  two  additional  types  of  spec¬ 
tra  giving  a  total  of  seven  types.  The  three  parameter  choice  made  in  this 
paper  is  justified  by  noting  that  Newton's  law  of  motion  relates  the  force  act¬ 
ing  on  a  vehicle  to  the  acceleration  input.  Therefore  as  far  as  vehicle  dynam¬ 
ics  is  concerned  there  are  three  stochastic  kinematic  quantities  to  be  obtained 
acceleration,  velocity  and  displacement.  This  suggests  a  three  parameter  sur¬ 
face  roughness  model  because  the  standard  deviations  of  velocity  a '  and  accel¬ 
eration  a  can  be  related  to  the  standard  deviations  of  the  slope  a  and  cur- 
vature  Oc,  respectively,  of  the  terrain  elevation  profile  (Part  III).  There¬ 
fore  the  calculation  of  the  basic  kinematic  quantities  associated  with  the 
vertical  dynamic  response  of  vehicles  (o^,  a v,  a )  requires  the  three  measured 
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roughness  quantities  0^,  og  and  0^  from  which  three  power  spectrum  parameters 
(C,  0  and  E)  can  be  determined. 

THREE  PARAMETER  MODEL  OF  DETRENDED  TERRAIN  ROUGHNESS  DATA.  This  section 
develops  the  formalism  and  mathematical  details  required  to  evaluate  the  power 
spectrum  of  a  measured  elevation  profile  from  the  RMS  values  of  the  elevation, 
slope  and  curvature  obtained  from  the  corresponding  detrended  elevation  pro¬ 
file.  This  is  done  for  each  of  the  five  types  of  power  spectra  described 
above . 


For  roughness  power  spectra  of  Types  1,  2  and  5  it  is  impossible  to  use 


values  of  o 


and  0^  obtained  directly  from  undetrended  data  to  calculate 


the  power  spectrum  coefficients  because  the  integral  P.dQ  diverges  (at  ft  -  0) 

o  d 


for  the  form  of  power  spectrum  given  in  equation  (12).  These  three  types  of 

terrain  roughness  must  be  handled  using  detrended  data  and  the  values  of  0^, 

0  and  o  for  detrended  data  can  be  related  to  integrals  of  the  form 
s 

j  P  <|>-2dft  which  converges.  For  roughness  spectra  of  Types  3  and  4  it  is 
o  d 


not  necessary  to  use  detrended  data  because  the  integrals  defining  o^,  0g 

and  0  are  defined  from  a  lower  cutoff  frequency  ft0  and  so  the  singularity  at 
c 

ft  =  0  is  avoided  and  undetrended  data  can  be  used.  But  the  detrended  formalism 
can  be  applied  to  spectra  of  Types  3  and  4,  and  this  section  evaluates  the  in¬ 
tegrals  that  appear  in  the  mathematical  expressions  for  0^,  0g  and  a c  (for  de¬ 
trended  data)  for  all  five  types  of  terrain  roughness  power  spectra. 


For  the  power  spectrum  given  by  equation  (12),  the  values  of  P^,  Pg  and 

P  for  undetrended  data  are 
c 


P,  =  Cft-2  +  Dft"3  +  Eft"4  (14) 

a 

Pg  =  (2n)zft2(Cfi'2  +  Dft'3  +  Eft'4)  05) 

Pc  =  (2n)4ft4(Cft*2  +  Dft'3  +  Eft“4)  [Sinn^AL"]4  06) 


The  form  for  the  slope  and  curvature  PSD  is  developed  in  Reference  1. 

The  necessary  integrals  for  the  detrended  scheme  of  the  three  parameter 
model  will  now  be  evaluated.  .This  will  be  done  for  each  of  the  five  basic 
types  of  roughness  power  spectrum. 

Type  1  Power  Spectrum 

The  variance  integrals  of  P.  and  Pg  cannot  be  performed  directly  on  (14) 
and  (15)  because  they  are  singular  for  Sft  =  0.  In  order  to  calculate  C,  D  and 
E,  the  variance  of  the  detrended  displacement,  slope  and  curvature  profiles 
must  be  used.  For  the  case  of  detrended  data  the  displacement,  slope  and 
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curvature  power  spectral  densities  are  P^0-2  >  Pg<>“2  an<*  Pc0“2  respectively 

where  $  is  defined  by  (7).  From  (14)  -  (16)  the  RMS  values  of  the  detrended 
terrain  displacement,  slope  and  curvature  are 


Jd  =  Vd  =  _-T0Dp^<t>~2dQ  =  Cfj(X)  +  Df2(A)  +  Ef3(A) 


a2  =  V  =  f°PJ-2  dfi  =  Cg,(A)  +  Dg2(A)  +  Eg3(A) 


a2  ~  V  =  /"p  dr2dfl  =  Ch^A)  +  Dh2(X)  +  Eh3(X) 

C  C  o  c 


The  values  of  ct^,  a2  and  a 2  are  represented  as  a  linear  combination 

of  integrals  where  the  f . (X) ,  g.(A)  and  h.(A)  are  given,  with  b  =  2n\  and 
a  =  rcAL,  as  1  1  1 


fi(X)  =  /  n-2(tr2da  = 


=  b  4/“— £ 

J  (1  + 


n2d» 

+  b2Q2)' 


f2(x)  =  /(TV2  dn 


00 

=  b4  f  — - 

J  (1  + 


b2fi2) 


f3(X)  =  j fi"4«))‘2dn  =  b 4J  - - 


b2ft2) 


gi  (X)  =  Un): 


00 

f  /sin  afl\ 

J  \  ,n  ) 


%-*dn  .  mpi 

a4 


/fl2  si 
(1  + 


in2 (afl)dfl 

+  b2n2)z 


oo  oo 

t,(A)  =  <2-o2  J  a-'  J  2-iia 


2(afl)dfl 

b2o2)2 
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g,(JD  =  <2«)*  j  !>-2  (^rfr2^  =  *^2^  J  siai 


(afi)dfi 

- z 

b2Q2) 


(25) 


MM 


=  «*)•  J  n2  J 


fl2  sin4  (afl)dfl 

(1  +  b2Q2)^ 


(26) 


h2  (A) 


-  «*>«  /  °  (^r9)4*-®  •  ^ 


(aO)dfl 

b2fi2)Z 


(27) 


h3(\)  =  (2n) 


00 

4  f  /sin  aflV 

J  \  aQ  / 


d>~2dQ  = 


(2n)4b4 

~7* — 


/sin4 
(1  + 


(aft)dQ 


(28) 


Note  the  integrals  gj,  g2 ,  hlf  h2  and  h3  diverge  where  a  =  0,  so  that  CTg  and 

ctc  diverge  for  AL  =  0.  Only  for  the  case  AL  f  0  are  o^,  ag  and  ctc  properly 

defined  by  infinite  integrals  when  0  -*  00  .  Also,  the  integrals  representing 
and  would  diverge  in  the  low  frequency  limit  if  the  detrending  factor 

<p~2  were  not  inserted  into  these  integrals.  In  other  words,  whereas  AL  f  0 
avoids  singularities  at  high  frequencies,  the  detrending  factor  <)>^2  avoids 
singularities  at  low  frequencies. 


The  values  of  the  integrals  in  (20)  -  (28)  are  evaluated  in  Reference  1 
by  complex  integration  using  the  contours  shown  in  Figure  2c,  and  have  the 
following  values 


MM  *  f 

(29) 

MM  =  jr 

(30) 

MM  =  ^ 

(31) 
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«,CM  ■  ^  [j(>  -  *  |  «-“/b 


(A)  =  [.»*  E,  (2a/b)  ♦  »‘2a/b  E*(2./b)j 

rw  n*b*Tl  (y  -2a/b\  a  -2a/b| 

*3(X)  =  Sr  [2  ('  ’  c  )  ‘  b  e  J 


.  rxx  _  4;i5b 

MM  ~  —4“ 


'3  .  r“/b  « 
8  8 


-2a/b  +  a  - 2a/b 
g  e 


(l  -  |  e-2*'b)] 


h2(M  =  |e2a/bEj(2a/b)  +  e‘ 2a/bE*(2a/b)  -  \  [e4a/bE, (4a/b) 

+  e‘4a/bE*(4a/b)]|  (36) 


MM  = 


4n5b3 


3  1  -2a/b/t  1  -2a/b\  a  -2a/b/  _  1 
8  ”  2  e  \X  ~  4  /  b  e  ^  2 


where  Ej  (x)  and  E*(x)  are  exponential  integral  functions  which  are  tabulated 
in  the  literature.  The  exponential  integrals  are  defined  as  follows 


Ej(x)  =  -Ei(-x) 


00 

'/v 


E*(x)  =  Ei(x) 


■If 


For  small  values  of  x  these  functions  have  the  following  series  expansions 
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E  (x) 

l 


v2  v3  v4  y5 

-Y  -  logeCx)  ♦  x  -  m  ♦  3T3  -  474  *  sTs  -  ••• 


X2  ,  X3  X4  .  X5 


E*(x)  =  Y  +  loge(x)  +  ^  +  2!2  +  3!3  +  474  +  5!5  +  **’ 


where  y  =  Euler  constant  =  0.5772157. 

Type  2  Roughness  Spectrum 

Detrended  data  must  be  used  for  this  case, 
and  are  (see  Figure  2b) 


The  integrals  defining  V^, 


“0 

°l=l 


P.<ir2dfi  =  CAfjCX.Oo)  +  DAf2(A,ft0)  +  EAf3(A,00) 


«o 

“•'I 


Psr2dn  =  CAg1(A,Q0)  +  DAg2(A,Q0)  +  EAg3(A,fi0) 


“0 

*  -f  V 


2dQ  =  CAhj (A,n0)  +  DAh2(A,Q0)  +  EAh3(A,fi0) 


where  the  integrals  Afj(A,ft0)>  Agj(A,Qo)  an<^  Ahj(A,fl0)  are  given  by 


Af^(A,fio)  = 


/fl2dfi  b  /n  .  n  n  \ 

- 2  =  x  10o  -  sin  60  cos  0OJ 

(1  +  b2fi2)  2  V  ' 


Af2(A,n 0)  =  b4 


“0 

hi 


(1  +  b2f)2) 


=  -Z-  sin2  0O  'v 
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M3(X,fl0>  =  b< 


Wo 

I  -  m  -  T  (°°* 

J  (1  +  b2ft2)  2  ' 


sin  6n  cos  B, 


o)^  b4Qo 


where 


00  =  tan"1 (b00) 


Agj (A,Q0) 


■  (~Jb4f  if; 


sin2(a0)dfl 
+  b2fi2) 


where 


=  |(^)2b4  [lo(Oo)  -  l[23)  (n0)  +  ii2a)  (fi0) 


L0(fi0)  =  fig  (5  -  §  b2ng  +  |  b4ng  -  |  b60g  +  ^  b«ng  -  ••  J 


L(2a)  ^  ^  _  flg  sin  (2afl0)  +  20n  cos  (2afln)  _  2  sin  (2aQc 

23  (2a)2  ‘  (2a)3 


li2a)  (fi0) 


=  2b2  <  4(2a)zftg  - 


J  flp  cos  (2afl0) 
J  (2a)4 


+  f(2a)4fig  -  12(2a)2ng  +  24I  -*ip  -^?*nC 
L  J  f2a) 


Ag2(\,fi0)  = 


/2tt\2i4  f  n  sin2  (aft)dfi 

w  J  „  ,  bW)‘ 

0 

=  \  (r)2fe4  [wo(Ao)  -  w[2a)(n0)  +  w^24)(00) 
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■ 


s 


where 


W0(00)  =  (i  ‘  b20g  +  b40g  -  b6ftg  +  b80g  -  •••) 


u(2a)  ,n  .  On  sin  (2aOn)  .  cos  (2a00)  1 

W1  (°<>)  2a  ( 2a f2a 


(2a)3  (2a)i 


W^2a)  (00)  =  2b2 


JfeoL- 

(L(2aF 


Ag3(^,flo) 


(2^F  (2a) 


■  (?)v  f  ^ 


(ifpjco8  (2afio> +  (2ir*  +[ll '  ufp] fio  sin  (2aQo)j 


2  (aO)dO 

- T~ 


(51) 


b202) 


=  |(f?)2b4  [*><*>)  -  «l [2a)(fto)  +  42a)(00)  -  •••] 


where 


H0 (00)  =  Oq^I  -  §  b20g  +  |  b40g  -  y  b60g  +  |  b80g  - 


H(2a)  (fio)  .  sin_i2aM 


H^2a)  (Q0) 


» 

„9  [og  sin  (2a0o)  200  cos  (2aQ0)  _ 

=  2b  L  2a  +  (2a)2 


2  sin  (2aflo) 
(2a)3 
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Ah  i  (A  ,Q(i )  = 


where 

.  (<a) 


44a)  C«0)  =  2b- 


Ah2(^,Wo)  = 


where 

WC«) 


fl2  sin4  (afl)dfl 
(1  +  b2fi2)* 


l  (F)4b4  |l  Lo(Qo)  “  2  [L»2a)  Wo)  -  b22a)(«o)  +  “•] 
+  \  ^4a)(n0)  -  44a)(fi0)  +  ii4a)(fi0) 


(52) 


/•«  \  _  sin  (4afio)  *  2fin  cos  (4aft0)  2  sin  (4afi0) 

tJ/°j  "  2a  +  (4 a)2  ‘  (4a  )a 


ft 


4(,4a)2;}(5 


«,"]  %  cos  (4aQo) 

2  J  (4.)* - 


♦  [(4a)««J  -  12l4a)20S  ♦  24^  S1°JfeB||)j 


/2n\4.  4  I  Q  sin4  (aO)dO 

\a  /  J  (1  +  b2n2)* 


(53) 


i  (?)v  {! 


W0(fi0)  -  2 


[«S 


2a 


'  (Qo) 


W*8a)(n0) 


(4a) 


(4a). 


(fl0)  "  V*&o)  + 


,)  = 


JTq  sin  (4aflp) 


4a 


cos  (4afl0)  _ ! 

(4a)2  (4 U2 
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w|4a)(fl0)  =  2b2  -  ^fyr]  cos  (4a00)  +  ^  Q0  sin  (4aft0)j 


Ah 3(\,fi0)  = 


_  /27t\4L4  f  sin4  (afi)dfi 

_wb  J  ;t;tw 

=  |  (r)4b4||  Ho  Wo)  -  2  ^H[2a)(n0)  -  42a)(Oo)  +  •••] 


\  [h^43)  C00)  -  44a)(Q0)  +  “•  | 


where 


„(“>«!„)  =  U°-I^k2 


„(4a),^  ^  _  „v2  [*($  sin  (4afi0)  .  2ft0  cos  (4aft0)  2  sin  (4aQ0)l 

h2  (q0)  -  2b  ^ -  +  -  (4D-2 — - (4IF~^J 


Type  3  Roughness  Spectrum 

Either  the  detrended  or  undetrended  procedure  can  be  used  for  this  case. 
The  integrals  defining  V^,  Vg  and  for  the  detrended  case  are  (see  Figure  2b) 


a 


P.<r2dfi  =  Cf  j  (A,fl0)  +  Df  2  (\,fl0)  +  Ef3(\,fto) 


a 

1 


Ps0’2dn  =  Cg^A.^o)  +  Dg2(A,ft0)  +  Eg3(A,fi0) 


Pc<rzdfi  =  Chi(A,fl0)  +  Dh2(A,Q0)  +  Eh3(A,Q0) 
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b 


The  integrals  f.(A,120),  g  (A,ft0)  an<*  h  (A,Q0)  are  given  by 
J  J  J 


00 

/O-J-j 

— ~ 

0  +  b 


dn 


- 2  -  f .  (A)  Af.(A,Uo) 

b2ft2)  J  J 


g .  (.A,120) 


00 

=  (™)\<  f  ,t.M- 

\a  /  J  (l  +  b2Q2)  1 


Agj (A.Qq) 


h.(A,i20J  =  I  “ 


'a  /  J  (1  +  b2Q2)2  j 

flo 


r —  =  h.(A; 


Zih^(A,fi0) 


1 58) 


(59) 


(60) 


Type  4  Poh  ;r  spectrum 

Either  the  detrended  or  undetrended  scheme  can  be  used  tor  this  case. 
The  integrals  defining  V.,  V  and  V  tor  this  case  are  (see  Figure  2b) 

Q  S  C 


flo 


<p'2d'J  =  Cty;A,Q0ffM)  +  Df2(A,fi0.ni)  +  EfJ(A,ft0,Ox)  (61) 


Pff'2dW  =  Cg^A.rto.tti)  +  DgalA.no.fix) 


+  Lg3(A,tt0,ftt)  (62' 


« 
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I 


Pc<T2dQ  =  Ch1(X,fi0,fi1)  +  Dh2(X,n0,fti)  +  Eh3(X,f20,«i) 


The  integrals  f.  (A^o,^) ,  g^  (A,fi0 >Qi).  and  hjCX.fto.Oi)  are  given  by 


fjCX.Qo.Qj)  =  b4 


r1 

/  n3~J 

J  a  +  b 


-^—2  =  Af  .(A.QO  -  Af  (X,n0) 
b2n2)  J  J 


8j  (X,flo»^i) 


■  ft)’-  fJvfb 


Sin2  (aO)dO  _  A  ^ 
- —  =  Ag  (X.Qj) 

b2fl2)  J 


Agj(X,ft0) 


hj (X,Q0,Q1) 


r  • 

_  /2n\4b4  I  fl3~J  sin4  (afl 

'  a'  J  (i  +  b2n2) 


=  Ah  (X,Oi) 

r>2\  J 


Ahj (X,fl0) 


Type  5B  Power  Spectrum 

The  detrended  scheme  must  be  used  for  this  case. 

V..  V  and  V  are  (see  Figure  2b) 
d  s  c 


The  integrals  defining 


+4.[ 


Pd4>-2dfi  +  /  Pd<t»-2dn 


=  C6f,(X,no,ni)  +  D6f2(X,Q0,Oi)  +  E6f3(X,n0.Oi) 


=  C6gl(A, «(,,«!)  +  D6g2(A,Q0»fli)  +  E6g3(\,n0»fii) 


(69) 


-  C6hj (A.Qj ,fj j )  +  D6h2(A,flo>^i)  +  E6h3(A,Qo,Qj ) 


The  integrals  fif ^ (A,fi0 .fii ) .  *8j (A.flo.^i) 


and  6hj  (A,ri0,Qi) 


are  given  by 


fif.CA.ilo.Q*) 


M.(A,n0j  -  =  £j0\)  +  M.(A,fl0) 


Af  (A.I^) 


(70) 


6g.(A,Q0,fii)  =  Agj(A,Q0)  +  g^A.fij) 


gj(A)  +  Agj (A,C0)  -  AgjCA.Qi) 


(71) 


fih  (A ,«„,«,)  =  Ah  (A,n0) 


+  h . (A.fij ) 


h.  (A)  +  Ahj  (A,n0)  -  Ah  (A,Ot) 


(72) 


Three  Parameter  Model  of  Undetrended 
Terrain  Roughness  Data 


This  section  considers  the  use  of  undetrended  nonstationary  elevation 
profiles  to  determine  the  surface  roughness  power  spectrum.  If  the  terrain 
roughness  power  spectrum  has  E  <  0  as  in  cases  Type  3  and  Type  4  of  Figure  2b 
there  is  a  lower  cutoff  frequency  fi0  for  the  power  spectrum  and  the  diffi¬ 
culties  associated  with  for  integrals  for  fl  ■*  0  do  not  arise.  For  this  case 
undetrended  roughness  data  can  be  used  to  determine  the  power  spectrum  coef¬ 
ficients  C,  D  and  E,  and  a  detrending  constant  \  does  not  enter  into  the 
calculations.  Types  1,  2  and  5  power  spectra  do  not  have  a  lower  cutoff  fre¬ 
quency  and  therefore  must  be  obtained  from  detrended  data. 

Type  3  Power  Spectrum 

For  a  Type  3  roughness  power  spectrum  that  exhibits  a  lower  frequency  cut¬ 
off,  the  integrals  expressing  0^,  o  and  for  undetrended  data  in  terms 
of  the  power  spectrum  coefficients  are  (see  Figure  2b) 


J2  =  J  (cq-2  +  dq-3  +  Efi'4)dfi  =  Cf!(n0) 


Df,(ftn)  +  Efa(ftn) 


(73) 


+  Dfi'3 


r<) 


+  Efi*4)  sin2  (afi)dfi 


(74) 


=  Cgi(Qo)  +  Dg2(fi0)  +  ^83(^0) 


+  Dfi'3 


sinn 


(afi)dfi 


(75) 


=  Chi (fi0)  +  Dh2(fi0)  +  Eh3(ft0) 


where  O,,  a  and  a  are  obtained  from  undetrended  data, 
d’  s  c 
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The  integrals  to  be  evaluated  are 


t,U20) 


dfi 

Q3 


.  v  ■  U3. 

=l  rf* 

L 

<Y> 

(l n\2  /  s^n2  (a(>)dO 

gjl^oJ  =  L,  J  I  q’Z 
■fio 


...  .  / 2rt\2  f  sin2  (aH)dn 

k?(o0)  -  ^  J  I  q3 


sin2  (aQ)dfl 

o* 


qr* 

/zn  1  f  Si::4  (afl)dQ 


M°o>  -  (f5)4  I 

MOo)  =  (2s)4  f 


(8: 


(8/ 


The  integrals  f,(Q0),  f2(«o),  f3(«o),  ga^o),  and  g3(fl0)  diverge  in  the  long 
wavelength  limit  Qq  -*  0. 

The  evaluation  of  these  integrals  is  done  in  Reference  1  and  yields  the 
following  results 


f,(fi0)  =  fto'1 
f2(flo)  =  ^2_ 

f3(fio)  =  (8 

gi(Oo)  =  ^r(n0-  J22a))  (8 

g2(fio)  =  ^(2^  ■  I323))  (8 

g3Cfio)  =  ^  (3^  ‘  ^2a)) 


(9 


^2 (fl0) 


-  /?n V  ('L  i 

"  \a  )  \16  0 


1  (2a)  1  T (4a)\ 

Q?  ‘  2  13  +  8  13  ) 


(92) 


h3(00)  = 


2n\4  /  1  1  I  (2a)  ,  1  , (4a)\ 

2  4  T  8  U  / 


a  /  V  8  0 


(93) 


where  the  T's  ire  given  by 


42a)  =  >  +  2asi  (2aQ0) 

“r> 


(94) 


C 


(2a)  _  -os 


20a 


-  +  2a2Ci(2a00) 

Oo 


14 


I 


(>J.l)  ..  cos 


0r 


30, 


+  4asi(4a00) 


^ia)  ^  cos  (4aO,J  _  2a_sin_14aOQ)  +  8a2Ci(4afto) 

3  205  00 


-  '! .J?.1-"- (-?9d2  -  2aj.-C°s  (2aOa)  _  4  a38i(2aQ0)  (96) 

JOy  ^O..  3 


(97) 


(98) 


(•-.a)  cos  (4a0n)  2a  sin  (4a0o)  8a2  cos  (4a00)  32  _3_:  r/.,r. 

=  -  30V - m - 30q  3 3  sl(4aS'° 


(99) 


and  rfhe/e  wine  integral  and  cosine  integral  functions  are  given  by 
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si(t) 


t 

■/ 


sin  x  .  7t 

-  dx  = 


t3  t5  t7 
-  -  +  t  -  - — •  +  - —  -  - —  +  •  •  • 
2  3!  3  5 !  5  7 !  7 


(100) 


/ 


Ci(t)  =  |  S25_*  dx 


t2  .  t4  t6 


Y  +  lo®efc  '  2!2  +  m  ‘  tl6  + 


(101) 


where  y  =  Euler  constant  =  0.5772157.  This  completes  the  analytical  evalua¬ 
tion  of  the  integrals  f^(Q0),  gj(fi0)  and  hj(fi0). 


Type  4  Power  Spectrum 

A  Type  4  roughness  power  spectrum  exhibits  a  spectral  window  having  lower 
and  upper  frequency  cutoffs  fio  and  respectively.  The  integrals  defining 
CT^,  and  are  written  as  (see  Figure  2b) 


o2  =  _f  (CQ-2  +  DA'3  +  Efi-4)dfi  =  Cfj(fi0,fii)  +  Df2(fi0,fii)  +  Ef3(fi0,fii)  (102) 
fto 


c 


W  i "" 


+  Dfi*3 


+  DJI-3 


+  Efi'4)  sin2  (afi)dfi 

=  0gj(fio>^l)  +  0g2  (fio  >^1  ) 


+  Efi'4)  sin4  (afi)dfi 

=  Ch^fiojfii)  +  Dh2(fio»Hi) 


+  Eg3(fi0,fii) 


+  Eh3(fio,fii) 


(103) 


(104) 


where  a.,  a  and  a  refer  to  undetrended  data,  and  where  the  integrals 
d  s  c 

fj(fi0,fii),  g j (fio ,fi j ) ,  and  h^(fio,fii)  must  be  evaluated  analytically. 

By  writing  the  integrals  in  equations  (102),  (103)  and  (104)  as 

ftj  OB  06 

/=/-/,  the  following  values  of  these  integrals  are  obtained, 
fio  fio  Hi 


fj (Q0,ni)  =  f.(fi0)  -  fj(fii) 


(105) 
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(106) 


gj(fto.fti)  =  gjt^o)  “  gj(^j) 
h.($W*i)  =  h^tfio)  -  h  (fl,)  (107) 


where  the  tunclrons  l.tlij,  g  (Vi,  and 
through  (9JJ.  1  1 


h^(V)  are  given  by  equations  (85) 


Calculations  of  the  Parameters 
ol  the  Roughness  Power  Spectnun 


inis  section  gives  the  procedure  lor  explicitly  calculating  the  power 
spectrum  coefficients  C,  0  and  £.  and  also  the  cutoff  frequencies  Q0  and  Hi 
which  occur  m  Type  1  -  Type  5  power  spectra.  The  values  of  the  integrals 
fj(A),  g.(A),  iiVA),  and  AK(A,Q0)>  Ag^(A,ilo)»  and  Ah^(A,Q0)  for  the  de- 

trer’ed  re  sc ;  g.  ((},>),  and  h^(flo)  for  the  Type  3  undetrended  case; 

and  <  .  ,  t(fio»^t)»  &nd  h.(£lo*^|)  l°r  The  Type  4  undetrended  case  have 

i  J  .1 


now  been  evaluated.  These  integrals  are  used  to  calculate  the  coefficient: 
T,  n  and  F  ( and  Q0  and  Z! ,  where  necessary)  that  describe  the  power  spectre? 
of  fhp  "N.  r»l  terr  .  i»:  e1 ovation  profile 


T'v  dot remiing  proceudre  can  be  applied  whether  or  not  there  is  a  lover 
cutoit  frequency  Ho*  arid  therefore  this  procedure  can  be  applied  it  <11  of  ';*ia 
five  pow«»*  sne<-*ruip.  types  shown  in  Figure  2b.  The  detrending  ptreteaure  must 
be  epyi’-ed  to  power  spectra  of  Types  1,  2  and  5  because  these  cases  have  pos- 
iti  ■  •  .slues:  of  power  spec  trust  in  the  limit  Q  *  0.  For  Types  3  and  4  power 
7n'  t,-i  cither  the  detrended  or  1  he  undetrended  formalism  can  be  used.  Firrt 
tb*  ..•Icilation  o."  C,  L  and  E  for  ail  five  types  of  terrain  roughness  >  i 1  I  tie 
dem  ur  iiv>  the  detrended  scheme,  and  then  C,  D  aid  E  vili  be  calculated  f  ir 
spe-.  .’a)  types  3  and  4  using  the  undetrended  scheme. 


*''»t  >  endvu  Sc  heme 


fvpu  1  (C  >  0,  D  >  0,  K  >  0) ,  Figure  2b. 


"V'e  tliree  simultaneous  equations  (17)  -  (19)  yield 


(V  f g (A)  f3( A) 
|V®  *Z(A)  g3(.\) 
h;I(A)  h3  (A) 

oft (A) . 


fl(A)  V  13(A) 
gi (A)  V°  g3(A) 
hi (A)  V®  h3(A) 

'detTa) 


(10?) 
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f 1 (A)  f2(A)  V. 

f 1 (A)  f 2 ( A)  f3(A) 

*i(A)  g2(M  v° 

ht (A)  h2(A)  V* 

DET(A)  = 

gi(A)  g2 (A)  g3(A) 
(A)  h2(A)  h3(A) 

DET(A) 


(109) 


which  give  C,  D  and  E  in  terms  of  the  measured  variances  V,,  V  and  V  of 
detrended  data.  S  C 

Type  2  (C  <  0,  D  <  0,  E  >  0) ,  Figure  2b. 

The  solution  of  the  three  simultaneous  equations  (17)  -  (19)  give 
the  power  spectrum  coefficients 


V  Af2(A,Q0)  Af3(A,fl0) 

V  Ag2(A,fto)  Ag3 (A jOo ) 
V*  Ah2(A,Q0)  Ah3(A,fi0) 

DETA(A.Qo) 


as 


D  = 


Af,(A,G0)  V  Af3(A,fi0) 
Agi(A,n0)  Va  Ag3(A>fi0) 
Ahi(A,fi0)  Vs  Ah3(A,Q0) 


DETA(A,Q0) 


(110) 


E  = 


Af  1  (A,fi0)  Af2(A M  V 
Agl(A,fi0)  Ag2(A,Q0)  V° 
Ah,(A,«0)  Ah2(A,ft0)  V* 

DETA(A,G0) 


DETA(A,Q0) 


Af ! (A,G0) 
Agi (A,fi0) 
Ah j (A,n0) 


Af2 (A,fi0) 
Ag2(A,Q0) 
Ah2 (A,Q0) 


Af3(A,fi0) 

Ag3(A,Q0) 

Ah3(A,ft0) 


'(HD 


and  the  equation  for  Q0  as 


Qq 


-D  -  >/D^~4CE 
2C 


(112) 


These  four  equations  determine  C,  D,  E  and  Q0. 

Type  3  (C  >  0,  D  >  0 ,  E  <  0) ,  Figure  2b. 

The  solution  of  the  three  simultaneous  equations  (17)  -  (19)  give  the 
following  expressions  for  the  power  spectrum  coefficients 
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V  f2(A,00)  f3CA,fio) 

V  82 (A  »do )  g3(A,ft0) 

Vc  h2(A,00)  h3(A,Q0) 


DET(A,Q0) 


fi(A.fio)  v  f3(x,n0) 

gl(A,fi0)  V°  g3(A,ft0) 

h,(A,«o)  h3(A,00) 

DET(A,Q0) 


(113) 


fi(A,«0)  f2(A,fi0)  va 
gi(A,Q0)  g2(A,fi0)  \ 
hj (A,fi0)  h2(A,0o)  Vc 

DET(A,»o)  ~ 


DET(A,»0) 


f  1  (A‘,fi0)  f2(A,«0)  ^(Aiflo) 
8i(A,00)  g2(A,00)  g3(A,Qo) 

hj(A,Qo)  h2(A,Q0)  h3(A,fl(j) 


(114) 


and  the  following  equation  for  Q0 


-D  +  J  D 


(115) 


These  tour  equations  determine  the  four  unknowns  C,  D,  E  and  fi0-  The  vari¬ 
ances  V,,  V  and  V  refer  to  data  detrended  with  the  constant  A. 
d’  s  c 

Type  4(C<0,  D>0,  E<0),  Figure  2b. 

The  three  simultaneous  equations  (17)  -  (19)  yield 


V  f2(A,«0,fi,)  f3(A,ft0,fi1) 

g2(A,n0,fl1)  g3(A,fl0.«i) 

Vc  h2 (A ,fl0 ,$1  j )  h3(A,fi0,fl1) 
D  ET  (  A,iT (7,q7)  ~ 


f i(A,Qo»^t)  V,  f3(A,n0.Oi) 

gi(A,fto>^i)  v  g3(A,fl0,ni) 

hitA.Qo.fli)  Vc  h3(A,flt>fli) 

DET(A,fi0,Q1) 


(116) 


f^A^o.Qt)  f2(A, fto.fit)  V 
g,(A,n0.(ii)  82(A,ft0,fi,)  V 
ht(At»o.Oi)  h2(A,n0,Q1)  Vc 

'  DEf(A7n^,oJ 


ii(A,llo)^j)  ^2(A,fto*(li)  f3(A,fic,fii) 
nuTfA  oot-  gi(A»00,n»)  g2(A,n0,Qi)  g3(A,Q0,ft1) 
u  ’  O’  W  h1(A,fi0,Q1)  h2(A,n0,01)  h3(A,Q0,n1) 


(117) 


u 


and  the  following  equations  for  the  bandwidth 


/D*  -  4CE 
2  C 


(118) 


-D  -  J D“ 


(119) 


These  five  equations  determine  the  five  unknown  power  spectrum  parameters  C, 
D,  E,  ft0  and  Qx. 

Type  5B  (C  >  0,  D  <  0,  E  >  0),  Figure  2b. 

The  three  simultaneous  equations  (17)  -  (19)  yield 


V  6f2(\,fi0,«i)  «f3(*,Oo.fli) 

V  6g2(A,f20,fi1)  6g3(X,Q0,fii) 
V®  6h2(X,n0»fli)  6h3(X,fl0.«i) 

DET(X^Wl7) 


6f 1(X,Q0,Q1)  V  6f3(X,ft0,Oi) 
6gi(x,o0,fli)  v°  6g3(x,fi0,n1) 

6hx(X,fi0>Q1)  VSq  6h3(X,Q0,«i) 
DET(X,n0>fii) 


(120) 


SfjCA.flo.fli)  6f2(A,Oo.«i)  vd 
68l(X,fio,fii)  6g2(A,«o,«i)  V 
6hl(X,fi0,fi1)  6h2(X,n0,«i) 

DET(X,fl0,fiiT 


6fi(x,n0,n1)  6f2(x,a0,«i)  Sf3(x,ft0,ni) 

nvTrx  no}-  ^8i(A»^o>^i)  6g2(X,Q0,Q|)  ^g3(A»fto*^i) 
DET(X,n0,Hjl  -  6hi(A>fio>fii)  6h2(Xrfi0,H1)  6h3(X,00,Q1) 


(121) 


The  equations  for  ft0  and  are 


■D  -  J Dz  -  ACE 


(122) 


-D  +  J Dz  -  4CE 


(123) 


These  five  equations  determine  C,  D,  E,  Qq  and  Qx . 
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For  the  detrended  scheme  the  values  of  a,,  a  and  o  are  obtained  from 

d  s  c 

data  measured  at  intervals  of  AL  and  detrended  with  a  factor  X.  As  such  a^, 
o  and  ct  represent  data  that  have  been  conditioned  in  order  that  the  theoreti¬ 
cal  determination  of  C,  D  and  E  can  be  accomplished  by  equations  (17)  -  (19). 
The  expressions  for  $“2  and  sin  (afl)/aQ  that  occur  in  the  right  hand  sides  of 
equations  (17)  -  (19)  account  for  the  detrending  factor  X  and  the  measurement 
interval  AL,  so  that  in  principle  the  coefficients  C,  D  and  E  are  independent 
of  the  intervals  of  data  measurement  or  the  detrending  technique,  and  reflect 
the  actual  roughness  condition  of  the  terrain.  It  is  clear  then  that,  as 
written  in  equation  (14)  through  (16),  describes  actual  surface  roughness, 
while  Pg  and  P^  still  contain  a  factor  containing  the  arbitrarily  selected 

value  of  AL.  But  AL  can  be  set  to  zero  at  this  point  of  the  calculation  to 

obtain  the  values  of  P  and  P  for  the  actual  terrain. 

s  c 

This  concludes  the  calculation  procedure  of  the  power  spectrum  coeffi¬ 
cients  C,  D  and  E  from  detrended  elevation  profile  data. 

Undetrended  Scheme 

Only  Type  3  and  Type  4  power  spectra  can  be  handled  using  undetrended 

data . 


Type  3  (C  >  0,  D  >  0,  E  <  0) ,  Figure  2b. 

For  a  Type  1  power  spectrum  the  undetrended  elevation  profile  can  be 
used,  and  the  power  spectrum  coefficients  are  determined  by 


V  f2(H0)  f3(«o) 
V,  g2(«o)  g3(fto) 

V®  h2(ft0)  MQ0) 

DET(ft0) 


fi(H0)  V.  f3(H0) 
gi(fio)  V°  g3(fl0) 
hi(fio)  V*  h3(fi0) 

DET(Q0) 


fl(ft0)  f2(fto)  V 
g,(fi0)  g2(Oo)  V° 
h^Qo)  h2(ft0)  \ 

DET(Oo) 


DET(Q0) 


fl(Ho)  f2(fi0)  f3(^o) 
gl(«o)  g2(Ho)  g3(Ho) 
hi (Ho)  h2(fio)  h3(Qo) 


(124) 


(125) 


and  the  following  equation  for  Q0 


Ho 


-D  +  VDZ  -  4CE 
2C 


(126) 


These  four  equations  determine  the  four  unknowns  C,  D,  E  and  fi0  that  describe 
a  Type  3  power  spectrum.  The  variances  V^,  Vg  and  Vc  refer  to  undetrended 
elevation  profile  data.  '  L 
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Type  4  (C  <  0,  D  >  0,  E  <  0),  Figure  2b. 

For  an  undetrended  elevation  profile  with  a  Type  4  power  spectrum  of 
bandwidth  flj  -  fl0,  the  power  spectrum  is  calculated  by 


V.  f 2  C^O  *dl  ) 

V  g2 (^0 )  g3(^0»^l) 

V*  h2(fi0,fli)  h3(Q0,n,) 
■  DET(fi0,fti) 


fi(Oo.Oi)  V.  f3(«o.Oi) 
g,(fto,fti)  V°  g3(Oo,Oi) 
h^Qo.ft,)  V*  h3(fl0,ni) 

DETOWh) 


(127) 


gi(Ho>°i)  82(»o,Qi)  Vs 
hi(Ho>^i)  h2(^o»^i) 

DET(ft0.«i) 


f2(Oo,Oi)  f3COo.fii) 

ni?Tro  o  X  -  8i(^o*^i)  82(fio.fii)  83(Oo.Oi) 

DET(00.«i)  -  hl(Q0.fii)  h2(00,nt)  h3(Q0,fti) 


and  the  following  equations  for  the  bandwidth 


/D*-  4CE 
2C 


-D  -  V Dz  -  4CE 


These  five  equations  determine  the  five  unknowns  C,  D,  E,  ft0  and  0t  that  de¬ 
scribe  the  Type  4  power  spectrum.  The  variances  Vrf,  Vg  and  Vc  refer  to  an 
undetrended  elevation  profile. 


Dominant  Frequencies  of  a 
Terrain  Elevation  Profile 


Roughness  power  spectrum  Types  3,  A  and  5  have  a  dominant  wavelength  for 
the  undetrended  elevation  profile  that  can  be  calculated  from  equation  (12) 
to  be 


■3D  ±  J 9DZ  -  32CE 


1 


E 


E 


(131 


For  detrended  elevation  profile  data  all  five  spectral  types  will  have  a  domi¬ 
nant  frequency.  An  approximate  expression  that  gives  the  dominant  wavelength 

of  detrended  elevation  data  in  terms  of  o.  and  a  is 

d  s 

=  k  * 2"  r  (132) 

D  s  • 


This  expression  is  affected  by  the  detrending  process. 


Numerical  Analysis  of  Terrain 
Elevation  Data 


A  numerical  study  was  done  to  determine  the  types  of  terrain  roughness 
power  spectra  that  are  associated  with  actual  terrain  areas,  and  to  determine 
their  basic  frequency  content.  Computer  program  RRFN  (Reference  1)  was  devel¬ 
oped  to  calculate  values  of  o^,  0g  and  a from  a  measured  elevation  profile. 

The  computer  program  RRFN  performs  the  following  specific  functions: 

a.  It  detrends  the  measured  elevation  profile  for  some  choice  of 
value  of  the  detrending  parameter. 

b.  It  calculates  the  standard  deviations  of  displacement,  slope 
and  curvature  for  the  detrended  and  undetrended  elevation 
profile. 

Computer  prgram  TERR  (Reference  1)  was  developed  to  calculate  the  rough¬ 
ness  power  spectra  from  the  values  of  o^,  Og  and  that  are  supplied  by  com¬ 
puter  program  RRFN.  The  computer  program  TERR  contains  all  the  integrals  and 
mathematical  functions  that  appear  in  Part  II.  Computer  program  TERR  performs 
the  following  specific  functions: 

a.  It  calculates  the  power  spectrum  parameters  C,  D  and  E  and  cutoff 
frequencies  Q0  and  fii  from  detrended  values  of  a^,  0g  and  CT  for 
spectral  Types  1  through  5,  and  from  undetrended  values  of  8^, 

Os  and  CT^  for  spectral  Types  3  and  4. 

b.  It  calculates  the  frequency  dependence  of  all  five  types  of 
terrain  roughness  power  spectra. 

The  elevation  profiles  of  100  terrain  sites  were  examined,  and  values  of 

0^,  <jg  and  0c  were  calculated  using  computer  program  RRFN  with  a  detrending 

constant  A  =  10  ft.  These  values  of  a,,  a  and  0  were  used  as  input  for  com- 

d '  s  c  r 

puter  program  TERR  to  calculate  the  values  of  the  power  spectrum  coefficients 
C,  D  and  E  and  the  associated  cutoff  frequencies  Qq  and  Q]  that  describe  the 
spectral  characteristics  of  the  actual  terrain.  The  results  appear  in  Table  2. 
The  values  of  C,  D  and  E  that  appear  in  this  table  were  calculated  using  the 
detrended  data  formalism. 
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The  results  in  Table  2  show  that  the  signs  of  the  coefficients  C,  D  and  E 
can  be  positive  or  negative  and  the  terrain  sites  can  be  separated  into  the 
five  spectral  classes  described  earlier  and  shown  in  Figure  2b.  The  frequency 
of  occurrence  is  shown  in  Figure  2d  which  indicates  that  spectral  Types  3,  4 
and  5  are  the  most  common  spectral  types. 

Figures  in  Reference  1  show  typical  elevation  profiles  that  correspond  to 
the  five  types  of  roughness  power  spectra.  These  figures  tend  to  support  the 
conclusions  that  can  be  drawn  from  Figure  2b,  namely:  that  Types  1  and  2  tend 
to  contain  mainly  long  wavelengths  (20-60  ft)  of  relatively  large  amplitude, 
Types  3  and  4  tend  to  have  a  roughly  sinusoidal  pattern  of  medium  sized  ampli¬ 
tudes  with  medium  sized  wavelengths  (5-10  ft),  while  Type  5  contains  a  high 
frequency  component  (wavelength  less  than  5  ft)  of  relatively  small  amplitude. 

The  frequency  of  occurrence  of  spectral  types  shown  in  Figure  2d  probably 
depends  on  the  type  of  terrain  selected  for  this  study  in  the  sense  that  all  of 
the  areas  selected  were  test  sites  used  for  the  operation  of  vehicles.  Other 
types  of  terrain,  such  as  may  be  encountered  in  macroroughness  studies  like 
the  elevations  measured  across  a  valley,  may  produce  a  somewhat  different 
occurrence  frequency  diagram.  The  lengths  of  the  elevation  profiles  used  for 
this  study  were  generally  300  ft.  The  selection  of  different  profile  lengths 
may  also  change  the  frequency  occurrence  diagram. 

It  was  shown  that  spectral  Types  3  and  4  exhibit  a  lower  cutoff  frequency 
Q0  and  it  is  therefore  possible  to  calculate  the  values  of  the  power  spectrum 
coefficients  C,  D  and  E  using  undetrended  elevation  profile  data.  Table  3 
gives  values  of  C,  D  and  E  calculated  for  Type  3  power  spectra  by  both  the 
detrended  and  the  undetrended  formalisms.  The  values  of  C,  D  and  E  pre¬ 
dicted  by  both  methods  are  essentially  equal,  as  they  should  be  since  they 
refer  to  the  natural  terrain  roughness  and  their  values  should  be  independent 
of  the  method  of  calculation. 


The  relationship  of  the  values  of  Og  and  with  the  values  of  0^  appears 

to  be  different  for  each  of  the  five  classes  of  terrain  roughness.  The  rela¬ 
tionship  is  shown  in  Figures  3a  and  3b  where  it  is  seen  that  there  is  some 
tendency  for  the  data  to  fall  into  distinct  groups.  Type  5  has  the  highest 

values  of  a  /a,  and  o  /a.  due  to  the  high  frequency  component  associated  with 
s  d  c  d 

this  spectral  type.  Types  1  and  2  have  the  lowest  values  of  these  ratios  be¬ 
cause  they  contain  dominant  low  frequency  components. 


Figures  3c  through  3e  give  the  dependence  of  the  values  of  0^,  0g  and  0^ 

on  the  choice  of  detrending  parameter  A.  The  values  of  decrease  rapidly 

with  A-*1,  while  the  values  of  o  and  0  decrease  more  slowly.  The  relatively 

s  c 

slow  variation  of  o  and  0  with  A-1  indicates  that  these  two  quantities  are 

associated  with  theShigher  frequency  components  of  the  terrain  profile,  while 

o.  is  associated  with  the  lower  frequencies, 
d 


The  values  of  the  spectral 
roughness  power  spectra  for  the 
given  by  equations  (14)  -  (16). 
are  shown  in  Figures  4a  through 


coefficients  C,  D  and  E  immediately  give  the 
terrain  displacement,  slope  and  curvature  as 
The  basic  forms  of  P^fl),  Pg(^)  an<*  PC(Q) 

6e  for  spectral  Types  1  through  5. 
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PART  III:  GEOMETRY  AND  KINEMATICS  OF  THE  VEHICLE  - 
GROUND  CONTACT 


Introductory  Remarks 


a 

In  the  previous  parts  of  this  paper,  terrain  roughness  was  described  in  a 
manner  independent  of  the  geometrical  characteristics  of  an  operating  vehicle. 
But  the  degree  of  roughness  of  a  terrain  area,  as  measured  by  the  dynamic  re¬ 
sponse  of  an  operating  vehicle,  will  depend  on  the  geometry  of  the  vehicle- 
ground  contact  areas.  Therefore  the  terrain  roughness  description  introduced 
earlier  in  this  paper  must  be  generalized  to  include  the  effects  of  the  geom¬ 
etry  of  the  ground  contact  area. 

Two  kinds  of  quantities  can  be  distinguished  in  vehicle  dynamics  problems' 
input  and  output  variables.  By  input  variables  are  meant  all  quantities  mea¬ 
sured  or  defined  at  the  soil-vehicle  contact  points.  Output  quantities  are 
all  quantities  measured  or  defined  at  points  on  the  vehicle.  A  theoretical 
prediction  of  the  output  dynamical  response  of  a  vehicle  requires  an  under¬ 
standing  of  interaction  of  the  geometry  of  the  vehicle  with  the  geometry  of 
the  ground  surface  because  this  interaction  will  produce  a  terrain  roughness 
power  spectrum  that  includes  the  effects  of  the  terrain-vehicle  contact  geom¬ 
etry.  This  terrain-vehicle  power  spectrum  serves  as  the  input  for  calcula¬ 
tions  of  the  dynamical  response  of  a  vehicle.  This  part  of  the  paper  calcu¬ 
lates  the  input  power  spectra  for  realistic  vehicle-ground  contact  geometries. 


Input  Power  Spectra  for  Single 
Point  Contact  with  the  Ground 


This  section  derives  the  power  spectra  associated  with  the  vertical  veloc 
ity  and  acceleration  of  a  point  which  is  constrained  to  move  along  a  specified 
elevation  profile  with  a  constant  horizontal  speed.  This  is  done  by  first 
calculating  the  standard  deviations  of  the  vertical  velocity  and  acceleration 
of  the  point  in  terms  of  the  standard  deviations  of  the  slope  and  curvature. 

The  standard  deviations  of  the  vertical  velocity  and  vertical  accelera¬ 
tion  of  a  point  that  follows  the  contours  of  the  surface  displacement  profile 
are  calculated  as  follows 
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where  Oy  =  standard  deviation  of  vertical  velocity  of  wheel-soil  point  of 
contact 

Og  =  standard  deviation  of  vertical  acceleration  of  wheel-soil  point  of 
contact 

v  =  vertical  velocity  of  wheel-soil  point  of  contact 
a  =  vertical  acceleration  of  wheel-soil  point  of  contact 

u  =  horizontal  velocity  of  vehicle  =  ^ 

At 

Equations  (133)  and  (134)  are  valid  for  both  detrended  and  undetrended 
data.  However,  they  are  valid  only  for  small  values  of  slope,  i.e.  for  the 
case  when  the  curvature  is  essentially  equal  to  the  second  derivative. 

Therefore  the  RMS  values  of  the  slope  and  curvature  immediately  determine 
the  RMS  values  of  vertical  velocity  and  vertical  acceleration  experienced  by  a 
point  travelling  over  an  elevation  profile  at  a  constant  horizontal  speed. 

This  is  done  by 


Ov  =  uos  (135) 

0a  =  u2Oc  (136) 

If  the  point  travelling  along  the  surface  elevation  has  an  associated  mass, 
the  standard  deviation  of  the  power  delivered  to  this  mass  is  approximately 
given  by 

a  =  ma  a  =  mu3a  o  (137) 

p  a  v  s  c 

where  a  =  standard  deviation  of  power 
P 

m  =  mass  of  point 

Although  the  values  of  a  and  o  are  of  some  interest,  they  cannot  be 
used  directly  to  calculate  the  corresponding  RMS  values  of  output  velocity  and 
acceleration.  The  calculation  of  these  output  quantities  requires  the  input 
power  spectra. 

The  power  spectral  density  functions  for  input  displacement,  slope,  cur¬ 
vature,  vertical  velocity,  and  vertical  acceleration  for  a  point  contact  will 
not  be  calculated.  The  PSD  functions  are  expressed  in  terms  of  the  spatial 
frequency  Q  or  the  time  frequency  f  that  are  related  as  follows 

Q  =  J  (138) 
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where  =  power  spectrum  of  input  velocity  and  Pg  =  power  spectrum  of  input 

acceleration.  Equations  (138)  through  (143)  can  be  used  to  determine  the 
power  spectra  in  terms  of  the  time  frequency. 

From  (133),  (134)  and  (138)  through  (143)  it  follows  that 


Pv(Q)  =  u2  Ps  (Q)  (144) 

P  (0)  =  u4  P  (0)  (145) 

3  C 

Using  (138)  it  is  clear  that 

P^(f)  =  u"1  Pd(Q)  (146) 

P’(f)  =  o'1  P  (Q)  (147) 

s  s 

P^.(f)  =  u_1  Pc(Q)  (148) 
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(149) 


P'(f)  =  u"1  P  (ft)  =  uP  (ft)  =  u2  P'(f) 

V  V  s  s 

P'(f)  =  u'1  P  (ft)  =  u3P  (ft)  =  u4  P'(f)  (150) 

a  a  c  c 

A  summary  of  the  one  point  power  spectra  is  given  as  follows.  Let 


»Jj(ft)  =  Cft-2  +  Dft"3  +  Eft-4 

(151) 

* °©3  ♦  E(?)4 

(152) 

Pd(ft)  =  *(0) 

(153) 

P^(f)  =  u-1  =  u*1  Pd(ft) 

(154) 

Ps(ft)  =  (2nft)2  Pd(ft) 

(155) 

p;(f)  =  (¥)2  pi(f) 

056) 

P  (ft)  =  (27tft)4  P.(ft) 
c  d 

(157) 

pc<f)  *  (¥)4  Pd(f) 

(158) 

P  (ft)  =  (2/luft)2  P.(ft) 
v  a 

(159) 

r(f)  =  (2nf )2  P^(f) 

(160) 

P  (ft)  =  (2nuft)4  P.(ft) 

*  a  a 

(161) 

P’(f)  =  (2nf)4  P’(f) 
a  a 

(162) 

Table  1  gives  the  units  of  the  various  power  spectra. 
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Input  Power  Spectra  for  a  Two 
Point  Contact  with  the  Ground 


The  power  spectra  for  the  case  of  two  contact  points  with  the  ground  can 
be  obtained  by  using  the  spectral  window  functions  given  in  equations  (15)  and 
(16).  The  various  kinematic  and  terrain  roughness  power  spectral  are  then  given 


II 

(163) 

pi<fi  -  ■■"*(;)  =  pd<(,> 

(164) 

Ps(Q) 

=  (2nd)*  [siniatl)]2  rd(0) 

(165) 

p;(f)  = 

(¥)2  pd«> 

(166) 

pc(0) 

=  (2nfi)4  [— a-^-a-n~]4  Pd(fi) 

(167) 

P^(f)  = 

(^)4  ^M^Y  pi«> 

(168) 

Pv(0) 

=  (2nufi)2  [sin  ^af!)]2  Pd(Q) 

(169) 

p;(f)  = 

(2„f)2  iMj  P.(f) 

(170) 

pa(fi) 

=  (2nufi)4  [S---aQ~^ '] 4  pd(Q) 

(171) 

p'(f)  = 
a 

(2nf),  [MAO]<  P.(f) 

(172) 

distance 

between  the  two  contact  points. 

The  one  point 

power  spectra  given  in  equations  (153)  through  (162)  are  regained  by  taking 
a  =  0  in  equations  (163)  through  (172),  i.e.,  the  spectral  window  functions 
have  a  unit  value  for  the  case  of  one  point  of  contact. 
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Spectral  Window  Functions  for  Vehicles 


C 

1 

1 


1 


R 


! 

f 


t 


The  calculations  of  the  input  power  spectra  for  vehicles  requires  a  knowl¬ 
edge  of  the  spectral  window  functions  for  vehicles.  The  slope  and  curvature 
spectral  window  functions  were  introduced  in  the  roughness  power  spectrum  cal¬ 
culations  to  account  for  the  finite  interval  of  measurement  of  the  elevation 
profile.  The  same  situation  arises  for  a  vehicle  operating  on  a  terrain  site. 
When  the  absorbed  power  or  vertical  acceleration  is  measured  at  the  driver's 
seat,  the  surface  roughness  is  indirectly  being  measured.  But  now  the  sam¬ 
pling  of  the  elevation  profile  is  being  done  by  the  wheels  of  a  truck  or  the 
track  of  a  tank,  and  two  physical  situations  must  be  considered.  First,  the 
sampling  is  done  with  finite  contact  lengths  and  second,  there  may  be  several 
unequal  intervals  of  measurement  as  in  the  case  of  varius  distances  between 
the  wheels  on  one  side  of  a  truck  running  over  an  elevation  profile.  Both  of 
these  physical  situations  must  be  described  to  calculate  the  spectral  window 
functions  for  vehicles. 

The  spectral  window  functions  can  be  obtained  by  averaging  over  all  pos¬ 
sible  two-point  intervals  along  the  areas  of  contact  with  the  ground.  The 
two-point  contact  spectral  window  functions  for  the  slope  and  curvature  are 
given  by  equations  (15)  and  (16).  Single  and  multiple  contact  areas  must  be 
considered  for  the  averaging  process  using  equations  (15)  and  (16). 

Single  Contact  Area 

For  the  case  of  a  single  contact  area,  such  as  produced  by  one  wheel  or 
by  a  tank  track,  the  spectral  window  functions  are  calculated  by  integrating 
equations  (15)  and  (16)  over  all  two-point  contact  positions  as  shown  in 
Figure  7a  with  the  result  that 


asw 


(2) 


Tr  T 


1  c  L  C L  sin2  [nft(x  -  y)J 

-  t r  J  J  dx  dy 


0  0 


(173) 


asw 


(4)  = 


T  T 
L  L 

/  / 

0  0 


sin4  [ rcfl( x  -  y)  1 


[nfi(x  -  y)  ] 


%  dx  dy 


(174) 


where 

asw  =  slope  spectral  window  function  for  single  contact  area 
asw  ^  =  curvature  spectral  window  function  for  single  contact  area 
T^  =  length  of  contact  area  in  the  direction  of  vehicle  motion 
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These  spectral  window  functions  are  evaluated  in  Reference  1  and  have  the 
following  values  for  small  values  of  QT^ 

asw<2>  (Tl)  =  1  +  a2(nTL)2  +  a4 (QT^)4  +  a6(0TL)6  +  ...  (175) 

asw<“)  (Tl)  =  1  +  02(OTl)2  +  M^)4  +  Pe(«TL)6  +  •••  (176) 

where 


«2 

18 

(177) 

«4  = 

:  *  ^ 

675 

(178) 

a6  = 

n6 

8820 

(179) 

02 

tt2 

9 

(180) 

04 

=  + 

75 

(181) 

06  = 

17n6 

13230 

(182) 

These  expressions  are  valid  for  Q  <  1/T^)  and  are  therefore  low  frequency 
approximations . 

For  large  values  of  fiT^  (high  frequency)  these  spectral  window  functions 
have  the  following  asymptotic  values 

asw  <21  •v  (183) 

1  L 

asw  <4)  •v  (184) 
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The  spectral  window  functions  for  the  case  of  a  single  contact  area  are  shown 
in  Figure  7b. 


Multiple  Contact  Areas 

The  calculation  of  the  spectral  window  functions  for  a  multiple  contact 
area  geometry,  such  as  is  produced  by  the  wheels  of  a  truck,  is  accomplished 
by  integrating  equations  (15)  and  (16)  over  all  two-point  contact  positions 
as  shown  in  Figure  7c,  with  the  result  that 


ASwl21  =  orribm; 


2  / 

i » j  L. 

itS 


T  T 

N  L  JL 


M  L.+T.  L.+T. 

N  l  L  j  L 


V*  r  r  si°2  "  H-)l 

Iff  [«n(L  -  nJP  dCi  dni 


0  0 


sin2  [7tQ(£.  -  r| j )  ] 

P(t  -  nJP  d^i  dnj 


i  j 


(185) 


ASV/(4)  (N  +  N !  /  2  )  T? 


T  T 
N  L  L 

2  /  / 


0  0 


„  L .  +T  L.+T. 
N  i  L  j  L 

I  f  / 

i ,  j  L .  L . 

.  ’ .  i  j 

itj 


sin4  [ 7TQ(f,  -  ni)l 

Inna,  -  n-)l4  dii  dni 


sin4  [nft(4.  -  n,)] 

Inna.  -  Oj)]5*  d^i  dnj 


(186) 


where 

ASW^2)  =  slope  spectral  window  (unction  for  a  truck 
ASV/(4)  =  curvature  spectral  window  function  for  a  truck 
N  =  number  of  wheels  on  one  side  of  a  truck 
=  position  of  wheel  i 
Lj  =  position  of  wheel  j 
^  ,r|j  =  integration  variables 

In  equations  (185)  and  (186)  it  is  assumed  that  all  wheels  of  the  truck  are 
the  same  size.  The  truck  spectral  window  functions  are  normalized  by  the  fac¬ 
tor  (N  +  N!/2)T2  to  give  a  unit  value  if  the  truck  dimensions  were  shrunk  to 
zero  size  to  produce  the  case  of  a  one  point  contact. 

The  truck  spectral  window  functions  defined  by  equations  (185)  and  (186) 
are  evaluated  in  Reference  1  for  the  three  types  of  wheel  geometries  shown  in 
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Figure  7d.  For  this  evaluation  it  is  assumed  that  the  wheel-ground  contact 
length  is  equal  to  the  tire  radius,  where  =  tire  radius. 


a.  Four  wheel  truck 


asw(2) 


■M 

■*{ 


2asw  (2) (T^)  + 


ASW  (4 )  =  i  <  2asw  (4  )  dR)  + 


[sin  (nLi 2fl)  1  ( 

'nh12k  J  j 

f sin  (nLlgn)*l  4) 

L  nL12ft  J  j 


b.  Six  wheel  truck 


asw(2.  -I  j3  .»(*)(V  *  2  ♦ 

f~sin  (7tL23ft)1  l 

L  nL23fl  J  j 


AS«(.)  .  i  j3  .-(‘XT,,  *  [^“]  ^  4 


sin  (nL23fl)  \ 

7XL23fi  | 


c.  Eight  wheel  truck 


ASW^2)  =  77T  <4  asw(2)(TD)  + 


■M 


fsin  (nL12n)~l  r 

L  nLi2^  J  L 


2  2 
+  fsin  (nLi3fl)*l 

^Ll3n 


•  ^  Z  r  "|  Z  p  - 

sin  (nL14fl)  sin  (nLgafl)  +  sin  (nL2<ift) 

ttL14Q  nL23fi  7iL24n 


sin  (nL 


r22]') 


(187) 


(188) 


(189) 


(190) 


(191) 


*"“>  -  fe  j4  ..«<*>(!„)  *  [•*%»£*>]  '  * 

f sin  (7iL14fl)~l  fsin  (nL2afl)"l  fsin  (nL24Q)~l 

J  |_  nL2  3fi  J  |_  J 

|'siii_(7tLaAn)"[  ) 

L  "L34«  J  f 


(192) 


where  L„  =  distance  between  wheels  i  and  j.  Figure  7e  shows  typical  averaged 
slope  and  curvature  spectral  window  functions  for  a  four  wheel  truck. 


Equations  (187)  through  (192)  show  that  the  truck  spectral  window  func¬ 
tions  introduce  the  relevant  geometrical  characteristics  of  the  wheel-ground 
contact,  including  wheel  spacing  and  tire  radius.  Equations  (175)  through 
(184)  show  that  the  relevant  geometrical  quantity  for  track-laying  vehicles  is 
the  track  length  in  contact  with  the  ground.  Ml  these  equations  can  be  easily 
generalized  to  the  case  of  half-track  vehicles. 


Contact  Length  Filter 


A  vehicle  will  not  respond  to  the  high  frequency  components  of  a  terrain 
elevation  profile  because  the  finite  contact  length  of  a  tire  or  a  track  with 
the  ground  tends  to  filter  out  these  frequencies.  For  instance  a  track  of 
length  Tl  will  filter  out  all  wavelengths  shorter  than  T^/2,  while  on  a  smaller 

scale  a  tire  of  radius  TR  will  filter  out  wavelengths  shorter  than  TR/2.  The 
situation  is  shown  in  Figure  8a. 

A  low  pass  filtering  action  due  to  the  finite  contact  length  with  the 
ground  is  physically  plausible  and  necessary,  but  a  rigorous  theoretical  deri¬ 
vation  of  the  form  of  this  filter  has  not  been  developed.  In  this  report  it 
is  assumed  that  the  filter  is  of  the  exponential  type 


-F(fi,TL) 

e 


(193) 


where  several  forms  of  the  filter  function  F(fi,Tj )  can  be 
Two  forms  of  the  function  F(fi,Tj  )  were  investigated 


a  a  T. 

F(Q,Tl)  =  «  “  ^  - 

s  s  V 


selected . 
in  this  report 


F(il,TL)  =  p(51TL)2 


(194) 

(195) 
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where  =  Q"1  =  spatial  wavelength.  For  the  form  appearing  in  equation  (194) 
which  is  linear  in  T^,  it  was  found  that  a  spatial  frequency  factor  °c/°g 
had  to  be  inserted  to  account  for  the  dominant  frequencies  of  the  terrain. 

Only  in  this  way  was  the  parameter  a  found  to  have  a  constant  value  indepen¬ 
dent  of  terrain  site  and  track  length. 

The  values  of  a  and  p  were  determined  empirically  from  experimental  ab¬ 
sorbed  power  data  for  trucks  and  track  laying  vehicles  (see  Part  V)  and  are 
equal  or  =  118.8  in.  (for  detrended  values  of  a  and  a  ),  a  =  123.2  in.  for 

undetrended  values  of  a  and  a  )  and  8  =  1.4.  SThe  form  of  the  filter  func- 

s  c 

tion  given  by  equation  (194)  appears  to  describe  the  absorbed  power  data  some¬ 
what  better  than  does  the  form  in  equation  (195),  so  that  all  further  calcula¬ 
tions  in  this  report  were  done  using  the  filter  function  represented  by 
equation  (194).  Figure  8b  shows  a  typical  exponential  filter  expressed  in 
terms  of 

The  filter  functions  can  also  be  written  in  terms  of  the  time  frequency 
using  equation  (138)  as 


(196) 

(197) 


where  t  =  T^/u  =  time  delay  for  transit  of  the  vehicle-ground  contact  length. 

Figure  8c  shows  a  typical  exponential  filter,  given  by  equation  (196),  ex¬ 
pressed  in  terms  of  the  time  frequency. 

INPUT  POWER  SPECTRA  FOR  VEHICLES.  Taking  the  spectral  window  functions 
and  the  low  pass  contact  length  filter  in  account  and  treating  vehicles  as 
rigid  bodies  gives  the  following  expressions  for  the  input  roughness  power 
spectra  and  the  input  kinematic  power  spectra  for  wheeled  vehicles: 


Pdi<!))  =  o'Fpd<°> 

(198) 

a  -F  (2) 

P  .(fi)  =  (2 nfi)2e  rASW  P.(Q) 
si  a 

(199) 

-F  (4) 

P  .  (fi)  =  (27ifl)4e  ASW  P.(Q) 
ci  a 

(200) 

„  -F  <2> 

Pv.(fi)  =  (2nuft)2e  ASW  Pd(fi) 

(201) 
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(202) 


-F  (4) 

P  . (Q)  =  (2nufl)4e  fASW  P.(Q) 
ai  d 


where  P^,  Pg^,  ^ci*  an<^  ^a£  =  input  vehicle  power  spectra  for  displace¬ 
ment,  slope,  curvature,  vertical  velocity  and  vertical  acceleration  respec¬ 
tively.  These  equations  are  valid  for  track  vehicles  with  ASW  ^  and  ASW  ^ 
replaced  by  asw  ^  and  asw  ^  respectively. 

Equations  (198)  through  (202)  shot/  that  the  input  kinematic  power  spectra 
for  vehicles  are  related  to  the  terrain  roughness  power  spectrum  that  has  been 
adjusted  for  vehicle  geometry.  The  input  kinematic  power  spectra  would  de¬ 
scribe  the  motion  of  a  vehicle  if  it  were  totally  rigid.  The  actual  vehicle 
response  is  calculated  as  a  joint  effect  of  these  input  kinematic  power  spec¬ 
tra  and  a  transmission  function  which  accounts  for  the  impedance  of  the  ve¬ 
hicle.  The  following  Part  IV  develops  the  transmission  function. 


PART  IV:  DYNAMIC  MODELS  FOR  MOVING  VEHICLES 


Introductory  Remarks 


It  is  of  value  to  vehicle  design  engineers  to  be  able  to  predict  the  ver¬ 
tical  dynamic  displacement,  velocity  and  acceleration  at  a  point  on  a  moving 
vehicle  and  ultimately  to  predict  the  power  absorbed  by  the  driver.  To  do 
this,  the  dynamic  response  of  a  vehicle  to  rough  terrain  must  be  known.  The 
dynamic  response  of  a  vehicle  can  be  described  by  a  transmission  function. 14-16 

Because  vehicles  undergo  pitching  and  rolling  motions  as  well  as  vertical 
motion  (heave),  the  power  absorbed  by  the  driver  will  depend  critically  on  his 
location  relative  to  the  center  of  mass  of  the  vehicle.  It  is  important  to  be 
able  to  separate  and  determine  the  contributions  to  the  absorbed  power  due  to 
the  surface  roughness,  the  internal  dynamics  of  the  vehicle,  and  the  location 
of  the  driver.  The  driver  location  variable  appears  in  the  pitching  motion 
transmission  function. 

This  paper  uses  simple  linear  viscoelastic  spring  models  to  describe  the 
dynamic  response  of  a  vehicle.  It  is  known  that  linear  viscoelastic  models  do 
not  adequately  describe  the  dynamic  properties  of  vehicles  because  the  vehi¬ 
cles  contain  nonlinear  springs  and  dashpots,  and  Coulomb  damping  in  addition 
to  viscoelastic  damping.  The  linear  viscoelastic  models  introduced  here  serve 
as  a  simple  expedient  way  to  evaluate  the  roughness  power  spectrum  method. 

The  transmission  functions  developed  will  be  used  in  Part  V  to  calculate 
the  power  absorbed  by  the  driver  using  a  power  spectrum  method.  The  absorbed 
power  calculation  requires  the  output  acceleration  power  spectrum  at  the 
drivers  seat,  and  this  can  be  calculated  using  a  transmission  function. 


Vertical  Motion  Transmission  Functions 


Single  Mass 
Vertical  Motion 

The  output  power  spectrum  of  a  mechanical  or  structural  system  is  related 
to  the  input  power  spectrum  through  the  transmission  function  for  the  system, 
Figure  8d.  The  transmission  function  of  a  dynamical  system  is  generally  a 
complex  number  whose  imaginary  part  is  a  measure  of  the  damping.  For  a  linear 
system  the  output  power  spectrum  is  given  by 


p;(f)  =  |T|2  P!(f) 


(203) 


where  P' (f) 
o 

T 


output  power  spectrum 
transmission  function 
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|T|  =  magnitude  of  transmission  function 
Pl^(f)  =  input  power  spectrum 

Mechanical  transmission  functions  generally  come  in  two  forms,  one  relating 
output  displacement  to  input  force,  and  the  other  relating  output  displacement 
to  input  displacement. 

If  the  input  to  a  linear  spring  system  is  a  force  and  the  output  is  a 
displacement  we  have  .  . 


P' 

do 


(204) 


where 


power  spectrum  of  the  output  displacement 

transmission  function  relating  input  dynamic 
force  to  the  output  dynamic  displacement 


P'  =  power  spectrum  of  the  input  force 


The  transmission  function  for  this  case  is4 


\ 


T 


dF 


_  1/k _ 

i  -  (f/f “P~+"l2ff7F 

n  n 


(205) 


(206) 


(207) 


where  k  =  spring  constant 

f  =  ■Jk/m/2n  =  natural  frequency 
n  _ 

£  =  C^/(2Vkm)  =  damping  ratio 
m  =  effective  mass 
Cp  =  damping  constant 
f  =  frequency 

* 

For  application  to  vehicle  dynamics  a  more  useful  transmission  (unction 
is  that  which  relates  output  displacement  to  input  displacement.  ...e  output 
power  spectrum  of  displacement  is  then  given  by 


Pdo(f)  = 


dd 


Pd,(f) 


(208) 
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where  T  =  transmission  function  relating  the  input  dynamic  displacement 
to  the  output  dynamic  displacement 

.  =  power  spectrum  of  input  displacement 
The  transmission  function  for  this  case  is  given  by4 


T  =  — 
d.i  1 


1  +  iqf 

(17  f  )z  +  iaf 
n 


=  T 


dd 


dd 


1  +  a2 f2 _ 

(f/fn)2T^  +  a2f2 


♦dd  =  tan_1 


af 3/f2 

_ n 

1  -  (f/f  )2  +  a2f2 

n 


(209) 

(210) 


(211) 


where  a  =  2^/1^.  It  is  easy  to  show  that  the  velocity-velocity  transmission 

function  and  the  acceleration-acceleration  transmission  function  satisfy  the 
following  conditions 


(212) 


where  T  =  transmission  function  relating  input  velocity  to  output 
velocity 

T  =  transmission  function  relating  input  acceleration  to 
output  acceleration 

Therefore  only  the  displacement-displacement  transmission  function  needs  to 
be  considered  for  vehicle  dynamics  problems. 

An  important  property  of  the  transmission  function  is  its  value  for  zero 
frequency 


Vf  = 0)  =  1 


(213) 


This  corresponds  to  a  simple  translation  of  whole  dynamic  system,  so  that  the 
static  output  displacement  equals  the  static  input  displacement.  Note  also 
that  jT^j2  ♦  f*2  wh.  i.  jT^pj2  *  f  4  high  frequency  limit.  The  trans¬ 

mission  function  for  a  vehicle  running  over  rough  terrain  will  depend  on  the 
vehicle  speed  because  (  the  relation  f  =  uil.  In  particular,  the  value  of 
j  T^ j  j 2  for  the  single  .egree  of  freedom  model  is  obtained  from  equation  (210) 

to  be 
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2  - 


(214) 


_ 1  +  o2u2fl2 _ 

(i  -  u2o2/f2)  +  o2u2n2 


Note  that  for  u  =  0  the  transmission  function  has  the  value  of 


1. 


Double  Mass 
Vertical  Hotion 

Vehicles  generally  exhibit  two  resonance  peaks  -  one  associated  with  the 
vertical  motion  of  the  body  (M  Hz)  and  the  other  with  the  vertical  motion  of 
the  wheel  suspension  system  (M5  Hz).  A  displacement-displacement  transmis¬ 
sion  function  for  a  two-mass  system  is  therefore  required  to  accurately  de¬ 
scribe  the  dynamic  response  of  a  vehicle  moving  over  rough  terrain.  The  deri¬ 
vation  of  the  two-mass  transmission  function  is  given  in  Reference  1,  and  is 
as  follows  (Figure  9a) 


T  _  A  +  iB 
dd  "  E  +  iF 


(215) 


where 


dd 


2  =  (AE  ♦  BF)2  +  (BE  -  AF)2 
(E2  +  F2)* 


A  *  1  - 


Vv 


CB  Cw 

B  =  w  I  t—  +  — 

B 


V 


urM, 


E  =  1  - 


W 


urM, 


-  u>2 


Vw  +  h 

,  kBkW  kW, 


F  =  tu 


S/  S  2 
^  kB‘w 


mb(cb  *  CW)  +  mwcb 
kBkW  kBkW 


where 


to 


MB 


2nf 

damping  constant  .for  vehicle  body 
damping  constant  for  wheel  suspension  system 
mass  of  vehicle  body  =  W^/g 
mass  of  wheels  =  W^/g 


(216) 

(217) 

(218) 

(219) 

(220) 


kD  =  spring  constant  of  vehicle  body 

D 

=  spring  constant  of  wheel  suspension  system 

It  appears  that  track  laying  vehicles  can  be  adequately  described  by  a 
single  mass  transmission  function,  while  wheeled  vehicles  generally  require  a 
double  mass  transmission  function  to  describe  their  vertical  motion.  Typical 
vertical  motion  transmission  functions  for  tanks  &nd  trucks  are  shown  in 
Figure  9b.  The  dynamical  parameters  for  several  track  laying  vehicles  are 
given  in  Table  4,  and  for  several  trucks  in  Tables  5  and  6. 


Combined  Pitching  and  Vertical  Motion 
Transmission  Function 


Because  a  vehicle  operating  on  rough  terrain  has  a  pitching  motion,  the 
values  of  the  power  absorbed  by  the  driver  are  expected  to  depend  on  the  geo¬ 
metrical  location  of  the  driver  relative  to  the  center  mass.  The  pitching 
motion  of  a  vehicle  can  be  described  by  a  transmission  function.  This  section 
calculates  the  displacement-displacement  transmission  function  for  a  combined 
pitching  and  vertical  mode  of  motion.  The  model  used  for  this  calculation  is 
shown  in  Figure  9c,  and  consists  of  a  rigid  rod  representing  the  vehicle  body 
that  is  supported  by  two  damped  springs  at  each  end. 

The  equations  of  motion  for  this  model  are  the  following  coupled  linear 
differential  equations'4 


mx  +  kjX  i  +  CjXj  +  k2x2  +  C2x2  =  kin,  +  +  k2r|2  +  C2r|2 


(221) 


J0  -  L1[k1xj  +  CjXi]  +  L2[k2x2  +  C2x2]  =  -L^kiflj 


+  Cjm )  +  L2[k2n2  +  c2n2) 


(222) 


where 

m  =  mass  of  vehicle 

x  =  vertical  acceleration  of  center  of  mass 
kj  =  equivalent  spring  constant  of  front  support  of  vehicle 
xt  =  vertical  displacement  of  front  support  of  vehicle 
Cj  =  equivalent  damping  of  front  support  of  vehicle 
Xj  =  vertical  velocity  of  front  support  of  vehicle 
k2  =  equivalent  spring  constant  of  rear  support  of  vehicle 
x2  =  vertical  displacement  of  vehicle  body  directly  above  rear  wheels 
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C2  -  equivalent  damping  constant  of  rear  wheels  and  suspension 
x2  =  vertical  velocity  of  truck  body  directly  above  rear  wheels 
r)i  =  elevation  of  ground  surface  at  front  wheels 
Hi  =  rate  of  change  of  due  to  passing  vehicle 
n2  =  elevation  of  ground  surface  at  rear  wheels 
n2  =  time  rate  of  change  of  r|2  due  to  moving  vehicle 
J  =  moment  of  inertia  about  pitch  axis  through  center  of  mass 
Lj  =  distance  from  front  of  vehicle  to  center  of  mass 
L2  =  distance  from  rear  of  vehicle  to  center  of  mass 

The  solution  of  the  coupled  differential  equations  is  given  in  Refer¬ 
ence  1,  and  is  given  in  the  form  of  a  transmission  function.  The  displacement- 
displacement  transmission  function  is  given  for  a  point  in  front  of  the  center 
of  mass  and  for  a  point  behind  the  center  of  mass  as  follows 


«)  *  (t!) 

(223) 

2  2 

(c?)  *  fe!) 

(224) 

where 

Tp  =  displacement-displacement  transmission  function  for  a  point 
a  distance  z  in  front  of  the  center  of  mass 

Tg  =  displacement-displacement  transmission  function  for  a  point 
a  distance  z  behind  the  vehicle  center  of  mass 

z  =  distance  from  center  of  mass  at  which  the  motion  is  to  be 
calculated 

and  where 


Cl  =  R  -  zS  (225) 
C2  =  T  -  zlJ  (226) 
C?  =  R  +  zS  (227) 
C2  =  T  +  zU  (228) 


A  ‘>9 


I 


R  =  (ZiW,  +  Z2W2)/(wi  +  W2) 

(229) 

1 

- 

S  =  (X,Ui  +  x2u2)/(u*  +  U2) 

(230) 

1 

T  =  (Z2W,  -  Z,W2)/(wi  +  W2) 

(231) 

• 

U  =  (X2Ux  -  U2X,)/(ui  +  u2) 

(232) 

Zx  =  aA(d  -  Jui2)  -  uwBg  -  bpA  +  whpfi 

(233) 

Z2  =  o(g(d  -  Jw2)  +  waAg  -  bpB  -  u)hPA 

(234) 

i 

W,  =  ik  -  n*ju2 )  (d  -  Jw2)  -  ui2Cg  -  b2  +  u)2h2 

(235) 

W2  =  ujC(d  -  Ju)2)  +  u)g(k  -  mu>2)  -  2u)bh 

(236) 

! 

Xt  -  PAWi  -  PbW2  -  bZx  +  u>hZ2 

(237) 

Xz  =  PBW,  +  PAWZ  -  bZ2  -  a)hZj 

(238) 

! 

Uj  =  W^d  -  Jw2)  -  wgW2 

(239) 

« 

U2  =  UJgW,  +  W2(d  -  Ju)2) 

(240) 

aA  =  kj  +  k2  cos  (wL/u)  +  ujC2  sin  (ioL/u) 

(241) 

► 

a  =  ioC i  -  k2  sin  (wL/u)  +  i uC2  cos  (wL/u) 

B 

(242) 

1, 

[ 

k 

\ 
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1 


eA 

pb 


L2[k2  cos  (u»L/u)  +  u>C2  sin  (tuL/u)]  -  Likt  (2A3) 

L2 [u)C2  cos  (u)L/u)  -  k2  sin  (u>L/u)]  -  ioCiLj  (244) 

L  =  Lj  +  L2  <245) 

b  =  k2L2  -  kiLj  (246) 

d  =  k,L?  +  k2L§  (247) 

h  =  (2L2  -  C,L,  (248) 

g  =  CjLf  +  C2L|  (249) 

C  =  Ct  +  c2  (250) 

k  =  ki  +  k2  (251) 


The  derived  transmission  function  can  be  used  to  calculate  the  motion  of 
a  vehicle  at  any  position  along  the  length  of  the  vehicle  relative  to  the  cen¬ 
ter  of  mass.  In  this  way  a  prediction  of  the  dependence  of  the  power  absorbed 
by  a  driver  on  his  location  relative  to  the  center  of  mass  can  be  made.  Fig¬ 
ure  9d  shows  a  typical  transmission  junction  for  the  combined  vertical  an 
pitching  modes  of  motion. 


PART  V:  ABSORBED  POWER  AND  VEHICLE  RESPONSE 


Introductory  Remarks 


Quantities  often  measured  for  a  vehicle  moving  over  rough  terrain  are  the 
power  absorbed  by  the  driver  and  the  vertical  acceleration  at  the  drivers  seat. 
The  power  absorbed  by  the  driver  is  the  energy  dissipated  as  heat.  It  repre¬ 
sents  the  energy  lost  due  to  vibration  damping  in  the  human  body.  Absorbed 
power  is  a  physiological  concept  that  has  been  developed  as  a  measure  of 
driver  fatigue  on  prolonged  exposure  to  vehicle  induced  vibrations.7  It  has 
been  found  empirically  that  the  maximum  power  that  a  driver  can  absorb  and 
still  function  reasonably  is  six  watts.  Therefore  it  is  of  importance  to  be 
able  to  estimate  the  vehicle  speed  at  which  six  watts  of  power  is  being  dissi¬ 
pated  by  the  driver.  The  six  watt  speed  clearly  depends  on  internal  vehicle 
characteristics,  the  location  of  the  driver  relative  to  the  center  of  mass  of 
the  vehicle,  and  the  surface  roughness. 

Two  methods  are  used  to  obtain  the  six  watt  speed.  The  first  is  a  formal 
power  spectrum  method  of  calculating  the  power  absorbed  by  the  driver  which 
includes  a  consideration  of  vehicle  dynamics.  The  second  is  a  regression 
analysis  of  the  six  watt  speed  measured  for  vehicles  operating  different  types 
of  terrain.  The  power  spectrum  method  introduces  the  terrain  roughness  and 
vehicle  geometry  through  the  power  spectra  described  in  Parts  II  and  III,  and 
the  vehicle  dynamics  through  the  transmission  functions  introduced  in  Part  IV. 
The  regression  analysis  uses  directly  the  detrended  values  of  ct^,  CTs  ai.d  Oc  to 

describe  the  terrain  roughness,  the  track  length  or  tire  radius  to  describe 
the  vehicle  geometry,  and  the  dimensionless  Froude  numbers  to  describe  the 
dynamical  characteristics  of  vehicles. 


Output  Power  Spectra  and  Vehicle  Response 


It  will  be  shown  subsequently  that  the  calculation  of  the  absorbed  power 
and  total  power  of  the  driver  requires  the  output  acceleration  and  velocity 
power  spectra.  The  power  spectra  of  the  vertical  motion  of  a  point  on  a  ve¬ 
hicle  can  be  obtained  from  equation  (208)  and  equations  (198)  through  (202)  to 
be 


P<!o(n)  *  ITJVV, 


(252) 


P  (Q)  =  (2nufl)2|  T .  .1  e'FASw(2)p  .(0) 
vo  I  dd|  d 

2  _r  . 

P  (0)  =  (2mift)4  It  .  .  e  Asw(4)p,(fi) 
ao  dd  d 


(253) 


(254) 


The  vehicle  speed  enters  directly  in  equations  (253)  and  (254),  but  the  de¬ 
pendence  on  the  vehicle  speed  is  also  indirectly  introduced  through  the 
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transmission  function,  as  for  instance  in  equation  (214).  The  corresponding 
power  spectra  expressed  in  terms  of  the  time  frequency  are  obtained  using 
equation  (154)  to  be 

Pdo(f)  *  |Tdd|2*'F!d<f>  (255) 

P’  (f)  =  (2nf)2  IT.  .1  2e'FASW<2>P'  (f)  (256) 

vo  I  ddj  d 

P'  (f)  =  (2nf)4|T.  .|2e"FASW<4>P' (f)  (257) 

ao  dd  d 


The  standard  deviations  of  the  vertical  displacement,  velocity,  accelera¬ 
tion  and  absorbed  power  at  a  point  on  the  vehicle  can  be  calculated  from  inte¬ 
grals  of  the  output  power  spectra,  given  by  equations  (252)  -  (254),  over  the 
frequency  ranges  of  the  five  basic  types  of  roughness  power  spectra.  Equa¬ 
tions  (252)  -  (254)  show  that  P^o  and  Pvq  diverge  for  f  =  0  (or  ft  =  0) .  This 

is  identically  the  same  problem  that  was  encountered  in  the  roughness  models 
of  Parts  I  and  II  where  it  was  shown  that  P^(ft)  diverge  for  ft  =  0.  Therefore 
the  same  difficulties  that  appear  in  the  roughness  models  of  spectral  types  1, 
2  and  5  appear  in  the  vehicle  response  problem.  Equations  (254)  and  (257) 
show  that  PaQ(ft)  and  PaQ(f)  are  not  divergent  for  ft  =  0,  corresponding  to  the 

situation  that  Pc(fl)  is  not  divergent  for  ft  =  0  as  described  in  Part  II. 

The  divergent  integrals  resulting  from  the  calculation  of  the  RMS  values 
of  the  output  vertical  displacement  and  velocity  can  be  evaluated  in  inserting 
a  filter  function  similar  to  that  defined  in  equation  (6).  Thus  the  following 
integration  factor  can  be  used 

0-2  =  UP^ft)!.  (258) 

[l  +  (27TY1ft)2J 


where  y1  =  filter  constant.  This  integration  factor  is  not  necessary  for  the 
calculation  of  absorbed  power  because  this  calculation  involves  only  the  ac¬ 
celeration  power  spectrum. 

The  standard  deviations  of  the  \ehicle  displacement,  velocity  and  acceler 
ation  are  given  by 


cr 


vo 


/ 
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»di)(n)*-2dn  =  / 


(259) 


(260) 
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(261) 


< 


h 


/OD 

P  (fi)dfl  =  /  P’  (f)  df 
ao  *&  ao 

The  integrals  in  equation  (259)  through  (261)  are  evaluated  only  in  the  re¬ 
gions  for  which  the  power  spectra  are  positive. 

4 

Absorbed  Power  and  Six  Watt  Speed  Calculated 
by  Power  Spectrum  Method 

The  absorbed  power  is  related  to  the  vertical  acceleration  of  the  driver 
(at  the  drivers  seat)  and  is  defined  in  terms  of  the  vertical  acceleration 
power  spectrum  as  follows6 


A  =  /  H2(f)P'  (f )df  =  J 

P  *tf  a  ao  *6 


A'  (f ,u)df 
P 


(262) 


where 

Ap  =  absorbed  power  (watts) 

Ha(f)  =  human  factor  function  jVwatts/ (ft/sec2 )j 

P'  (f)  =  output  vertical  acceleration  power  spectrum  at  the  location  of 
the  drivers  seat 

A^(f,u)  =  absorbed  power  spectral  function 
The  human  factor  function  H  (f)  is  empirically  determined  and  has  the  form  of 

3 

a  band  pass  filter.6  This  function  appears  in  Figure  9e. 

The  frequency  dependence  of  the  absorbed  power  is  given  by  equa¬ 
tions  (254),  (257)  and  (262)  as 


A  (Q,u)  =  (2nufl)4  IT.  ,|2e'FH2ASW(4>P.(0) 
p  |  da |  a  a 

A’(f,u)  =  (2rtf)4|T.  .]2e"FH2ASW(4>P'(f) 
p  |  ad |  a  d 


(263) 

(264) 


where  A^(Q,u)  =  uAp(f,u).  Equations  (262)  and  (264)  show  that  five  basic 

quantities  enter  into  the  absorbed  power  calculation,  and  these  depend  on  the 
vehicle  dynamics,  vehicle  geometry,  driver  location  and  terrain  roughness: 

a.  transmission  function 

b.  low  pass  filter  associated  with  the  ground  contact  length 

c.  human  factor  function 
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d.  curvature  spectral  window  function  for  vehicle 

e.  power  spectrum  of  terrain  elevation 


The  roughness  power  spectrum  is  given  by  the  three  parameter  model,  and  refers 
to  actual  terrain  described  by  the  parameters  C,  D  and  E  of  Part  II. 

Because  the  functions  and  H  that  appear  in  equations  (263)  and  (264) 

-F 

are  intrinsically  functions  of  the  time  frequency  f,  while  the  functions  e 
and  ASW^4)  are  intrinsically  functions  of  the  spatial  frequency  Q,  it  follows 
from  f  =  uft  that  the  absorbed  power  given  by  equation  (262)  always  depends  on 
the  vehicle  speed.  The  six  watt  vehicle  speed  is  obtained  from  the  absorbed 
power  equation  (262)  by  taking  A^  =  6  watts. 

It  should  be  pointed  out  that  the  absorbed  power  dissipated  by  the  driver 
is  not  the  same  as  the  total  power  associated  with  the  kinetic  and  potential 
energies  of  the  motion  of  the  driver.  The  RMS  value  of  the  total  power  de¬ 
livered  to  the  driver  is  given  by 


CT 

po 


P  (ft)Pa„(fl)dfi 
vo  ao 


(265) 


where 


ct  =  RMS  value  of  total  power  of  driver 
po 

Wp  =  weight  of  driver 

A  comparison  of  the  experimental  values  of  absorbed  power  with  values 
predicted  by  the  numerical  integration  of  equations  (262)  and  (264)  are  shown 
in  Figures  10a  through  11c  for  a  series  of  track  laying  vehicles.  Figure  lid 
gives  a  comparison  between  predicted  and  measured  values  of  the  six  watt  speed 
for  track  laying  vehicles.  A  similar  procedure  for  trucks  appears  in  Fig¬ 
ures  12a  through  12d. 


Linear  Regression  Predit tion  of  the  Six  Watt  Speed 


The  six  watt  speed  is  the  vehicle  speed  at  which  the  driver  is  absorb¬ 
ing  six  watts  of  power;  it  is  a  measure  of  ride  quality.  A  high  value  for 
the  six  watt  speed  means  a  smooth  ride,  and  a  low  value  means  a  rough  ride. 

The  value  of  the  six  watt  speed  depends  on  both  terrain  and  vehicle 
characteristics. 

For  the  purpose  of  the  design  and  testing  of  vehicles,  it  is  important  to 
have  a  simple  expression  for  the  six  watt  speed  giving  its  dependence  on  ter¬ 
rain  roughness  parameters  and  on  vehicle  parameters.  In  this  way  the  results 
of  vehicle  tests  conducted  at  different  terrain  areas  with  different  types  of 
vehicles  can  be  compared,  and  an  evaluation  of  design  changes  can  be  made. 

The  vehicle  parameters  describe  the  geometry  and  internal  dynamics,  while  the 
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terrain  parameters  describe  the  displacement,  slope  and  curvature  of  an  eleva¬ 
tion  profile. 

Values  of  the  six  watt  speed  for  several  vehicles  have  been  measured  for 
a  number  of  terrain  sites  where  elevation  profiles  have  been  determined.  The 
regression  analysis  given  here  represents  the  six  watt  speed  as  simple  powers 
of  the  terrain  roughness  parameters  and  the  vehicle  parameters.  This  analysis 
is  valid  only  if  the  six  watt  speed  is  a  monotonical ly  increasing  or  decreas¬ 
ing  function  of  the  chosen  parameters.  It  was  found  that  the  vehicle  dynamics 
parameters  k,  m  and  C  could  not  be  successfully  brought  into  the  regression 
analysis,  and  therefore  only  terrain  roughness  and  vehicle  geometry  were 
considered. 

A  regression  analysis  incorporating  k,  m  and  C  could  not  be  accomplished 
because  the  dependence  of  the  six  watt  speed  of  these  parameters  could  not  be 
represented  by  a  simple  power  law.  The  dependence  on  the  vehicle  weight,  for 
example,  was  found  to  be  either  direct  or  inverse  depending  on  the  choice  of 
the  resonance  frequency  or  the  damping  ratio  as  a  basic  vehicle  dynamics 
parameter;  both  choices  gave  essentially  the  same  quality  of  fit  to  the  mea¬ 
sured  six  watt  speeds.  This  suggests  that  the  dependence  of  the  six  watt 
speed  on  the  vehicle  weight  is  more  complicated  than  a  simple  power  law.  In 
fact,  the  six  watt  speed  is  expected  to  initially  increase  with  the  vehicle 
weight  producing  a  smoother  ride,  and  then  decrease  with  additional  weight 
giving  a  rougher  ride  due  to  hitting  the  bump  stops. 

It  is  possible  to  find  a  direct  correlation  between  the  six  watt  speed, 
the  three  terrain  roughness  descriptors  a,,  o  and  CT  for  detrended  data,  and 

Q  S  C 

the  vehicles  parameters.  Both  wheeled  vehicles  and  track  laying  vehicles  were 
considered  for  this  study. 

Track  Laying  Vehicles 

Seven  mathematical  forms  were  chosen  for  the  regression  analysis  of  the 
six  watt  speed  for  track  laying  vehicles.  They  are  as  follows: 


u6  =  A o“ 

(266) 

u6  =  Aa^ 

°  s 

(267) 

u6  =  Aa* 

(268) 

“S  =  TL 

(269) 

u„  =  Aa“  j[ 

(270) 
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u6  =  Act"  tJ  jJ  (271) 

“«  =  Jj  (272) 


where 

u6  =  six  watt  speed  of  vehicle 

°d’  V  °c  =  values  of  displacement,  slope  and  curvature,  respec¬ 
tively,  of  a  detrended  elevation  provile  (A.  =  10  ft) 

=  track  length 

=  jounce  length 

The  coefficients  in  equations  (266)  through  (272)  were  obtained  by  fit¬ 
ting  these  equations  to  measured  six  watt  speeds.  Eight  track  laying  vehicles 
were  used  in  this  analysis,  LEO  2AV,  M60  Al,  ATR,  SMI,  AISV,  HIMAG-5,  MICV  and 
the  Ml 13.  The  six  watt  speed  was  measured  for  these  vehicles  for  a  total  of 
twenty-four  terrain  sites.  The  dimensions  of  the  quantities  shown  in  equa¬ 
tions  (266)  through  (272)  are  chosen  to  be  as  follows:  [u6]  =  mph,  [a^]  =  in, 
[as]  =  1,  [acJ  =  in.*1,  [Tj]  =  in.  and  [Jj  =  in.  * 

This  data  combined  with  a  library  regression  computer  program  gave  the 
following  results  for  the  six  watt  speed  and  the  coefficient  of  fit  (cf): 


u6 


cf 


“6 


23.55 

d 

0.29 

9.62 
CTo  rim 


cf  =  0.02 


(266a) 


(267a) 


-  *7. 2 

un  -  Oo . z n 

c 

cf  =  0  04 


(268a) 


0.393  T°'874 

uf>  = - -0.875 - 

d 

cf  =  0.70 


(269a) 
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(270a) 


3.33  J°98 


cf  =  0.82 


1.324  TO  -289  jO  .761, 

L  L 


(271a) 


0.944  a0706  T?-412  J°-672 
s  L  L 


(272a) 


Equations  (266a)  through  (268a)  show  that  the  detrended  displacement  o. 
is  the  relevant  terrain  roughness  parameter  for  track  laying  vehicles,  while 
Og  and  ac  play  secondary  roles.  Equations  (269a)  through  (271a)  show  that  the 
jounce  length  is  the  primary  vehicle  parameter  for  track  laying  vehicles 
while  the  track  length  plays  a  secondary  role.  The  best  fit  to  the  mea¬ 
sured  si*  watt  speed  data  is  given  by  equation  (272a)  with  all  parameters  in¬ 
cluded,  but  the  fit  is  not  much  better  than  given  by  equation  (270a)  with  only 
the  two  parameters  and  J^.  A  comparison  between  the  experimental  six  watt 
speed  and  those  predicted  by  equation  (272a)  is  shown  in  Figure  13a. 

When  the  vehicle  dynamics  parameters  k,  m  and  C  are  eventually  brought 
into  the  expression  for  u6,  they  should  be  entered  as  nondimensional  parameters 
in  the  form  of  the  damping  ratio  and  the  Froude  numbers.17 

Wheeled  Vehicles 

Seven  mathematical  forms  were  chosen  for  the  analysis  of  the  six  watt 
speed  for  wheeled  vehicles.  They  are  expressed  in  terms  of  the  tire  radius 
(T^)  and  wheel  base  length  (Lyg)  as  follows: 


*  o 

U6  =  Aad 


(273) 


u6  =  Aa* 


(274) 


u6  =  A o' 


(275) 


.a  6 

U«  =  AOs  LWB 


(276) 
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Ug  =  Act 

D  s 

Te 

R 

(277) 

a 

u6  =  Act 
°  s 

T£ 
wB  JR 

(278) 

u6  =  A od 

ctP  aY  lA  Tp 
s  c  WB  R 

(279) 

The  coefficients  in  equations  (273')  through  (279)  were  obtained  by  fit¬ 
ting  these  equations  to  measured  six  watt  speeds  for  trucks.  Five  trucks  were 
used  in  this  analysis:  PAC-CAR,  TARADCOM-HMTT,  10  ton  (8x8)  -  Cargo  truck  A; 
DRAGON  WAGON,  10  ton  (8x8)  -  Cargo  truck  B;  GERMAN  MAN,  10  ton  (8x8)  -  Cargo 
truck  C;  M656,  5  ton  (8x8);  M520E1  GOER,  8  ton  (4x4);  and  the  Czeckoslovakian 
TATRA  813. 


A  library  regression  routine  (MULFIT)  was  used  to  analyze  the  six  watt 
speed  data  for  these  five  trucks,  and  gave  the  following  results  for  the  six 
watt  speed  and  the  coefficient  of  fit  (cf): 


14.87 

u6  -  JTT.3T? 


cf  =  0.38 


0.041  1 
U6  =  (jl~7T2 


cf  =  0.51 


0.789 

u6  -  jrm 


cf  =  0.62 


1285 


cf  =  0.76 


(273a) 


(274a) 


(275a) 


(276a) 


217.9 


cf  =  0.80 


(277a) 


317.7 


cf  =  0.80 


(278a) 
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(279a) 


1099  a°-475  a0942 
a  c 

u6  “  n2TZ2  tCT.  ISB  T J  .  59 
as  Kra  tr 

cf  =  0.82 

Equations  (273a)  through  (275a)  show  that  the  detrended  slope  o  is  the 

•  ® 

significant  terrain  roughness  parameter  for  wheeled  vehicles,  while  o  and  CT^ 
play  subordinate  roles.  The  best  fit  to  the  measured  six  watt  speeds  of 
wheeled  vehicles  is  given  by  equation  (279a)  with  all  variables  included,  but 
this  fit  is  not  much  better  than  shown  in  equation  (277a)  with  only  two  param¬ 
eters  a  and  T^.  A  comparison  between  the  experimental  values  of  the  six  watt 
speed  and  those  predicted  by  equation  (279a)  is  given  in  Figure  13b.  The 
relatively  poor  agreement  shown  in  Figure  13b  suggests  that  an  alternative  ve¬ 
hicle  parameter,  such  as  the  jounce  length,  is  required  to  obtain  a  better 
correlation  between  measured  and  predicted  six  watt  speeds  for  wheeled 
vehicles. 

Dimensional  Analysis 

Dimensional  analysis  uses  relationships  between  nondimens ional  quanti¬ 
ties  which  remain  valid  irrespective  of  the  scale  of  the  physical  system.17 
The  procedure  is  to  calculate  a  set  of  relevant  dimensionless  parameters  for  a 
physical  system.  For  this  study  the  physical  system  is  a  vibrating  vehicle 
moving  over  rough  terrain. 

The  Froude  numbers  for  a  physical  system  are  given  by17 


(280) 


where  £  =  Froude  number,  M  =  characteristic  mass,  V  =  characteristic  speed, 
F  =  characteristic  forces  operating  in  the  system,  and  L  =  characteristic 
length. 

The  Froude  numbers  associated  with  a  linear  harmonic  oscillator  are 


,  MV2  _  MV 
"  CVL  CL 


(281) 


(282) 


where  f^  =  resonance  frequency  of  vehicle,  K  =  spring  constant  of  vehicle,  and 
C  -  damping  constant  of  vehicle.  Because  the  resonance  frequency  of  a  vehicle 
is  a  directly  measurable  quantity,  equation  (281)  is  used  as  the  appropriate 
Froude  number  and  therefore  the  appropriate  dimensionless  parameter  for  repre¬ 
senting  the  six-watt  speed  for  vehicles  is 
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(283) 


t  =  _“fi_ 
U*  fRTL 


For  wheeled  vehicles  the  characteristic  lengths  would  be  the  tire  radius  and 
wheelbase  length.  The  other  dimensionless  parameter  for  a  vibrating  system  is 
the  damping  ratio 


(284) 


A  third  parameter  that  was  examined  is  the  ratio  of  jounce  length  to  ground 
contact  length 


k  =  jl/tl  (285) 

For  a  vibration  problem  it  is  the  vehicle  mass  M  =  W/g  that  enters  (281) 
through  (284)  not  the  vehicle  weight.  Dimensionless  parameters  such  as  W / (KL) 
and  W/(CV)  are  inappropriate  for  a  vibration  problem.  The  correct  dimension¬ 
less  parameters  equivalent  to  (283)  are  MV/(CL)  and  CV/(KL).  This  can  be  seen 
from  the  extreme  case  of  zero  gravity  where  g  =  0  and  W  =  0  but  M  =  some  fixed 
value;  in  this  case  vibrations  can  still  occur.  On  the  other  hand  for  a 
statics  problem,  such  as  deformation  of  a  solid  (vehicle)  under  the  action  of 
gravity,  it  is  the  weight  and  not  the  mass  that  enters  the  appropriate  dimen¬ 
sionless  parameters. 

Terrain  roughness  has  been  described  by  the  three  parameters  ct^,  ag 
and  o  .  Within  this  description  the  terrain  roughness  dimensionless  param¬ 
eters  are 


l 


s 


=  o 

s 


(286) 


a  o , 
c  d 


(287) 


A  combination  of  these  parameters  that  was  examined  is 


^R 


c  d 


(288) 


A  number  of  dimensionless  relationships  were  developed  with  corresponding 
coefficients  of  fit  from  data  for  eight  track  laying  vehicles,  and  are  as 
follows 


i  =  0.021  n  28  (289) 

u6  K 

cf  =  0.60 
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(290) 


(291) 


(292) 


Equation  (292)  shows  an  insignificant  dependence  on  the  ratio  of  jounce 
length  to  track  length.  A  comparison  of  measured  values  of  Ue/(fpTj)  with 
those  values  if  shown  predicted  by  equation  (291)  is  shown  in  Figure  13c. 
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PART  VI:  CONCLUSIONS 


CONCLUSIONS .  This  paper  generalizes  the  method  of  characterizing  the 
ground  surface  displacement  profile  by  a  one  parameter  power  spectrum  of  the 
form  Cft*2.  A  three  parameter  model  of  the  power  spectrum  of  the  terrain  dis¬ 
placement  profile  data  is  introduced  in  such  a  manner  that  the  additional 
parameters  will  more  accurately  describe  the  contribution  of  the  long  wave¬ 
lengths  to  the  surface  roughness,  and  will  describe  the  cases  where  periodici¬ 
ties  are  present  in  the  terrain  such  as  in  the  case  of  plowed  fields. 

The  three  roughness  power  spectrum  parameters  are  used  to  calculate  the 
power  spectra  and  standard  deviation  values  of  the  vertical  displacement,  ve¬ 
locity,  acceleration  and  power  absorbed  by  the  driver  for  a  vehicle  operating 
on  rough  terrain.  A  comparison  of  theoretical  and  experimental  results  was 
made . 


The  theoretical  and  experimental  studies  of  terrain  roughness  and  dynamic 
vehicle  response  yields  the  following  conclusions: 

a.  A  three  parameter  power  spectrum  model  can  be  used  to  classify 
terrain  roughness  into  five  basic  types  which  adequately  de¬ 
scribe  natural  and  manmade  terrain  features  including  periodici¬ 
ties  (Parts  I  and  II). 

b.  The  three  parameters  of  the  roughness  power  spectrum  can  be  de¬ 
termined  from  the  standard  deviations  of  the  displacement,  slope 
and  curvature  of  a  detrended  elevation  profile,  but  in  some 
cases  undetrended  data  can  be  used  (Part  II). 

c.  The  power  absorbed  by  the  driver  of  a  vehicle  depends  on: 

(1)  the  human  factor  function  relating  absorbed  power  to  the 
acceleration  power  spectrum,  (2)  the  vehicle  transmission  func¬ 
tion  that  describes  the  internal  dynamics  of  a  vehicle,  (3)  the 
low  pass  ground  contact  length  filter  that  describes  the  filter¬ 
ing  effect  of  a  track  or  wheel,  (4)  the  vehicle  spectral  window 
functions  that  introduce  the  geometry  of  the  vehicle-ground  con¬ 
tact,  and  (5)  the  three  parameter  roughness  power  spectrum  that 
introduces  surface  roughness  (Parts  III,  IV  and  V). 

d.  The  six  watt  speed  for  track  laying  vehicles  can  be  related  by  a 
regression  analysis  to  the  vehicle  geometry  and  the  three  rough¬ 
ness  descriptors  a,,  o  and  a  for  detrended  data.  The  signifi- 

d  s  c 

cant  terrain  roughness  parameter  for  track  laying  vehicles  is 

a.,  while  for  wheeled  vehicles  it  is  o  (Part  V). 
d*  s 
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Table  1 

Units  of  the  Power  Spectra 


c 

in. 

D 

dimensionless 

E 

in.  -1 

Pd(fl) 

in.  3 

P’(f) 

in.2  sec 

P  (fi) 

in.3  sec'2 

V 

p;(f) 

in.2  sec-1 

Pa(Q) 

in.3  sec’4 

p;(£) 

in.2  sec’3 

Ps(0) 

in. 

p;(f) 

sec  _1 

Pc(f>) 

in.'1 

P’(f) 

in.'2  sec 

c 
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Table  4 

Dynamical  Parameters  of  Track  Vehicles 


Vehicle 

Type 

Weight 

lb 

Spring  Constant 
lb/in. 

Damping  Constant 
lb  sec/in. 

Track  Length 

in. 

M113A1 

24,200 

6,060 

3110 

105 

MICV 

43,000 

8,400 

720 

150 

AISV 

2,940 

300 

80 

50.25 

HIMAG-5 

79,910 

10,000 

400 

168 

LEO  2AV 

111,000 

15,750 

8000 

188.3 

M60  A1 

106,000 

20,000 

8000 

171 

I 

I 


520 


Table  5 

Dynamical  Parameters  for  Trucks 


Wheel  Assembly 


Vehicle  Type 

Weight 

lb 

Body 
Spring 
lb/ in. 

Damping 
lb  sec/in. 

Weight 

lb 

Spring 
Constant 
lb/ in. 

Damping 
Constant 
lb  sec/in. 

PAC-CAR 
TARADCOM-HMTT 
10  ton  (8x8) 

46,179 

4718 

830 

5131 

75,490 

830 

DRAGON  WAGON 

10  ton  (8x8) 

44,136 

4509 

822 

4904 

72,150 

822 

GERMAN  MAN 

10  ton  (8x8) 

47,610 

4864 

887 

5290 

77,829 

887 

M656 

5  ton  (8x8) 

23,040 

2354 

250 

2560 

37,664 

250 

M520E1-GOER 

8  ton  (4x4) 

41,400 

2500 

1000 

•  - 

_  _ 

•  _ 

Table  6 

Geometrical  Properties  of  Trucks 


Tire 

Radius 

L,2 

1*13 

L14 

L23 

L24 

1*34 

Vehicle  Type 

in. 

in. 

in. 

in. 

in. 

in. 

in. 

PAC-CAR 

TARADCOM-HMTT 

10  ton  (8x8) 

26 

58 

190 

248 

132 

190 

58 

DRAGON  WAGON 

10  ton  (8x8) 

26 

58 

202 

260 

144 

202 

58 

GERMAN  MAN 

10  ton  (8x8) 

23.75 

65.5 

209.5 

275 

144 

209.5 

65.5 

M656 

5  ton  (8x8) 

22 

56.25 

147.75 

204 

91.5 

147.75 

56.25 

M520E1-GOER 

8  ton  (4x4) 

34.5 

235 

-- 

-- 

-- 

-- 

-  “ 

Pd(il)  (m2/cycle/m) 


n(CYCLES/METER) 
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Figure  11 
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